t  * 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


gbSTAT^ 


AD-A241  071 


DTIC 

,  -LECTE 
'  OCT  0  3  1991 


THESIS 


APPROXIMATE  INTERVAL  ESTIMATES  FOR 
MECHANICAL  RELIABILITY 


by 

Yang,  Wen-Huei 


September,  1990 


Thesis  Advisor: 


W.  Max  Woods 


Approved  for  public  release;  distribution  is  unlimited 


91-12189 

iiiniiiiif 


91  10  2  057 


Unclassified 


security  classifies''’  i  ol  this  page 


la  Report  Securiiy  Classification  Unclassified 


Ha  Security  Classification  Authority 


2b  Declassification  Downgrading  Schedule 


J  Performing  Organization  Report  Numbcr(s) 


t>a  Name  of  Performing  Organization 

Naval  Postgraduate  School 


be  Address  (dry,  stare,  and  ZIP  code ) 

Montcrcv.  CA  93943-5000 


8  a  Name  of  Funding  Sponsoring  Organization 


Sc  Address  (city,  state,  and  ZIP  code)' 


REPORT  DOCUMENTATION  PAGE 


lb  Restrictive  Markings 


6b  Office  Symbol 
(If applicable)  OR 


Sb  Office  Symbol 
(If  applicable) 


3  Distribution  Availability  of  Report 

Approved  for  public  release;  distribution  is  unlimited. 


5  Monitoring  Organization  Report  Number(s)  ' _ 


7  a  N  ame  of. Monitoring  Organization 

Naval  Postgraduate  School 


7b  Address  (city,  state,  and  ZIP  code) 

Monterey.  CA  93943-5000 


9  Procurement  Instrument  Identification  Number 


10  Source  of  Funding  Numbers  ! 


Program  F.Iemcnt  No  Project  No  Task  No  Work  Unit  Accession  No 


it  Title  (Include  scanty  classification  t  APPROXIMATE  INTERVAL  ESTIMATES  FOR  MECHANICAL  RELIABILITY 


12  Personal  Author(s)  Wcn-Huci  Yang 


13a  Type  of  Report 

Master's  Thesis 


13b  Time  Covered 
From  To 


14  Date  of  Report  fycar.  month,  day)  15  Page  Count 

September  1990  64 


16  Supplementary  Notation  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  policy,  or  po¬ 
sition  of  the  Department  of  Defense  or  the  U.S.  Government. 


17  Cosf.ti  Cocos 


Field 

Groun 

IS  Subject  Terms  I  continue  on  reverse  If  necessary  and  Identify  by  bloc 1'  number) 

Reliability,  Confidence  Limit,  Parametric,  Normal,  Unknown  Means  and  Variances 


19  Abstract  t  continue  on  reverse  if  necessary  and  identify  by  block  number; 

Two  approximate  interval  estimation  procedures  for  mechanical  component  reliability,  P(X>  I),  are  developed  and  their 
accuracy  evaluated  by  computer  simulations.  The  strength,  X,  of  the  component  and  the  stress,  Y,  applied  to  it  arc  inde¬ 
pendent  normally  distributed  variables  with  unknown  means  and  variances.  In  the  first  interval  procedure  the  variances  arc 
equal.  In  the  second  procedure  the  variances  may  be  unequal. 

The  derived  intervals  arc  quite  accurate  for  the  cases  simulated  which  include  large  and  small  sample  sizes.  These  proce¬ 
dures  are  simple  to  apply  and  require  the  use  of  percentile  points  of  the  Student's  t  distribution.  In  the  second  procedure,  the 
degrees  of  freedom  of  the  associated  t  statistic  is  a  function  of  the  test  data  ,  and  therefore  it  is  random. 


20  Distribution  Availability  of  Abstract 
E  unclassified  unlimited  □  same  as  report 


■22a  Name  of  Responsible  Individual 
W.M.  Woods 


DD  FORM  1473.S4  MAR 


21  Abstract  Security  Classification 
□  dtic  users  Unclassified 


22b  t  elephone  < include  Area  code) 

(408)  646-2768 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  arc  obsolete 


22c  Office  Symbol 

OR  Wo, 


security  classification  of  this  page 


Unclassified 


i 


Approved  for  public  release;  distribution  is  unlimited. 


Approximate  Interval  Estimates  for  Mechanical  Reliability 

by 

Wen-Huei  Yang 
Commander,  Taiwan  Navy 
B.S.,  Chinese  Naval  Academy,  1976 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
■  September  1990 

^u^lor:  _ 'v  ^  _ 

Wen-Huei  Yang 

Approved  bv:  [/J  .  x  _ 

W.M.  Woods,  Thesis  Advisor 


=y,  Second  Reader 


Peter  Purdue,  Chairman, 
Department  of  Operations  Research 


11 


X 


ABSTRACT 

Two  approximate  interval  estimation  procedures  for  mechanical  component  reli¬ 
ability,  P(X>  I-),  are  developed  and  their  accuracy  evaluated  by  computer  simulations. 
The  strength,  X,  of  the  component  and  the  stress,  Y,  applied  to  it  are  independent 
normally  distributed  variables  with  unknown  means  and  variances.  In  the  first  interval 
procedure  the  variances  are  equal.  In  the  second  procedure  the  variances  may  be  une¬ 
qual. 

The  derived  intervals  are  quite  accurate  for  the  cases  simulated  which  include  large 
and  small  sample  sizes.  These  procedures  are  simple  to  apply  and  require  the  use  of 
percentile  points  of  the  Student's  t  distribution.  In  the  second  procedure,  the  degrees 
of  freedom  of  the  associated  t  statistic  is  a  function  of  the  test  data  ,  and  therefore  it  is 
random. 
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I.  INTRODUCTION 


Let  X  and  Y  be  independent  random  variables  with  normal  cumulative  distribution 
functions  Fx(x)  and  Fr(y)  respectively.  Suppose  X  is  the  strength  of  a  mechanical  com¬ 
ponent,  and  Y  is  the  stress  applied  to  the  component.  The  strength  depends  on  the 
material  properties,  manufacturing  procedures  and  other  factors.  The  stress  is  a  function 
of  the  environment  to  which  the  component  is  subjected.  Component  failure  is  defined 
by  the  event  X  <  Y .  Component  reliability,  R,  is  defined  by  R  =  P( X  >  Y) .  R  is  called 
mechanical  reliability.  Two  lower  confidence  limit  procedures  for  R  are  developed  in  this 
thesis.  In  both  procedures  the  means  and  variances  of  X  and  Y  arc  unknown.  In  one 
procedure  the  variances  are  assumed  to  be  equal. 

A  nonparaimric  interval  estimation  procedure  for  R  was  first  proposed  by 
Birnbaum  [Ref.  1:  pp.  13-17)  using  the  Mann- Whitney  U  statistic.  Birnbaum  and 
McCarty  developed  a  procedure  for  computing  the  minimum  sample  size  needed  for 
such  a  confidence  interval  to  have  a  given  width  and  confidence  level  [Ref.  2: 
pp.  55S-562).  Owen,  Craswcll  and  Hanson  [Ref.  3:  pp.  906-924]  provided  more  detailed 
tables  for  use  in  computing  sample  sizes  and  confidence  intervals  for  the 
Birnbaum-McCarty  procedure.  Tables  designed  especially  for  the  normal  distribution 
were  also  included  in  their  paper.  Govindarajulu  [Ref.  4:  pp.  ,229-238)  observed  that  the 
bounds  employed  by  Birnbaum  and  McCarty  for  obtaining  the  confidence  intervals  can 
be  substantially  improved  asymptotically  and  reduced  the  Birbaum-McCarty  bounds  by 
approximately  1,2.  Church  and  Harris  [Ref.  5:  pp.  49-54]  pointed  out  that  the  sample 
sizes  required  by  these  nonparamctric  procedure  are  likely  to  be  too  large  for  many 
practical  situations. 


Gwen  [Ref.  6:  pp.445-47S]  gave  an  exact  confidence  limit  procedure  for 
R  =  P(X  >  x)  where  X  is  normally  distributed  with  unknown  mean  and  unknown  vari¬ 
ance,  and  .v  is  a  constant.  His  procedure  uses  the  noncentral  t  distribution  and  extensive 
table  lookups  arc  needed.  Owen  and  Hua  [Ref  7:  pp.  285-31 1J  developed  special  tables 
that  reduced  these  calculations.  Their  tables  were  limited  to  two  confidence  level  values; 
namely  90%  and  95%.  Lee  [Ref'S:  pp.  15-22]  reports  on  a  closed  form  equation  for  this 
confidence  limit  that  is  approximate  but  quite  accurate.  His  equation  applies  for  any 
confidence  level  and  uses  the  central  student's  t  distribution. 

Church  and  Harris  [Ref.  5:  pp.  49-54]  developed  approximate  confidence  limits  for 
R  =  P(X  >  Y)  ,  under  the  assumption  that  the  stress,  Y,  has  a  standard  normal  distrib¬ 
ution  and  Fx(x)  is  normally  distributed  with  unknown  mean  and  variance. 

Throughout  this  thesis  we  write  X  ~  S(n,  o-)  to  denote  that  X  has  a  normal  distrib¬ 
ution  with  mean  n  and  variance  g2. 


II.  APPROXIMATE  INTERVAL  ESTIMATION  PROCEDURE  FOR 
RELIABILITY  R  =  P(X  >  Y)  -  EQUAL  VARIANCE  CASE 

A.  LOWER  CONFIDENCE  LIMIT  PROCEDURE 

Suppose  the  strength,  X,  of  a  mechanical  device  and  the  stress.  Y,  applied  to  it  are 
independent  variables,  with  normal  probability  distributions.  We  assume  both  means 
are  unknown  and  both  variances  are  unknown  but  have  common  values  a2.  The  me¬ 
chanical  reliability,  R,  of  the  device  is  defined  as  follows: 


R  =  P[X>  >3 

X—  >'-  (i‘x  —  ,u  >•) 


=  P 


c\'~ 


<?v'2 


1‘X-J‘Y 
<rv2  . 


(2.1) 


where  <l>  is  the  standard  normal  cumulative  distribution  function.  Let  5—  a 
Then 


R  =  <i><; 


)- 


(2.2) 


A  consistent  estimator  {Ref.  9:  PP.2S9-294]  of  5  is 


K  = 


(2.3) 


where 
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7:  pp.  285-311]  have  taken  this  approach  for  the  univariate  case.  Their  method  required 
the  extensive  use  of  tables.  We  shall  take  a  different  approach  here,  in  order  to  find  a 

A 

closed  expression  for  an  approximation  to  dL. 

In  his  thesis,  Lee  [Ref.  8:  pp.  15-22]  developed  an  approximate  100(1  —  a)%  lower 
confidence  limit  for  R  =  P(X  >  x )  where  X  ~  iX(n,  a1),  n  and  a 2  are  unknown,  and  *  a 
constant.  Lee's  expression  is 


A 

It  1-3 


=  Kr 


Ki 


2  (n-yfc) 


m. 

fl-a,  n-1’ 


(2:5) 


where  Kx  =  ~  A  ,  and  r, is  the  100(1  -  a),h  percentile  point  of  the  student's  t  dis¬ 

tribution  with  n- 1  degrees  of  freedom. 

A  method  analogous  to  Lee's  procedure  can  be  used  to  develop  an  equation  for 

A 

in  lhe  bivariate  case  described  at  the  outset  of  this  section.  If  we  substitute 
(Hn  +  1/m) .for  (1  In)  and  (n  +  m  -  1)  for  (/i)  in  Equation  (2.5),  we  would  obtain  the  lower 
confidence  limit 


h,^  =  x- 


V  n  >»  ) 


K 


2(n  +  m  -  1  -  yjk  ) 


1/2 


0-a,  n-m- 2‘ 


(2.6) 


where  K=  .  -  .  The  corresponding  100(1  -<x)%  lower  confidence  limit  for 
R  =  P{X>Y)  is  : 


RL,  l-a 


(2.7) 


Our  computer  simulation  results  show  that  this  lower  confidence  limit  is  quite  accurate. 
The  results  are  tabulated  in  Tables  1  and  2.  A  description  of  the  simulation  procedure 
along  with  a  analysis  on  computer  results  will  be  given  later  in  this  chapter.  The  im- 

A 

portant  point  of  these  tables  is  the  comparison  between  R  and  Rm9[,.lhiM.,r  If  these  two 


5 


values  are  equal,  the  approximate  lower  confidence  limit  procedures  given  in  Equations 
(2.6)  and  (2.7)  are  nearly  exact. 
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Table  1.  ANALYSIS  OF  90%  CONFIDENCE  LIMIT  APPROXIMATION  OF 
R  =  P(X>Y)  FOR  EQUAL  VARIANCES  CASE  WITH  EQN.  2.6 


n 

m 

8 

8 

8 

30 

20 

30 

8 

8 

8 

30 

20 

30 

8 

8 

8 

30 

20 

30 

8 

8 

8 

30 

20 

30 

8 

8 

8 

30 

20 

30 

8 

8 

8 

30 

20 

30 

.8928 


.9047 


.9033 


.8928 


.9033 


.9450 


.9527 


.9515 


.9451 


.9527 


.9516 


.9S77 


.9908 


.9903 


.9S77 


.9908 


.9903 


.9130 


.8880 


.8890 


.9130 


.8880 


.8890 


.9130 


.8870 


.8900 


.9130 


.8870 


.8900 


.9200 


.S960 


.9200 


.8820 


.8960 


Mean  error 
from  R 

Variance  of 
error 

.1398 

.0110 

.0068 

.0704 

.0033 

.1398 

.0110 

.0990 

.0068 

.0704 

.0033 

.1092 

.0073 

.0722 

.0040 

.0508 

.0019 

.0073 

.0721 

.0040 

.0508 

.0019 

.0563 

jHEEEEHH 

.0310 

.0009 

.0212 

.0004 

.0563 

.0025 

.0310 

.0009 

.0212 

.0004 

Table  2.  ANALYSIS  OF  95%  CONFIDENCE  LIMIT  APPROXIMATION  OF 


R  =  P(X  >  Y)  FOR  EQUAL  VARIANCES  CASE  WITH  EQN.  2.6 


R 

n 

n 

m 

A 

■Rioooa-j),  id-i) 

True  confi¬ 
dence  level 

Mean  error 
from  R 

Variance  of 
error 

8 

8 

.8855 

.9650 

.1954 

.0126 

1.0 

8 

30 

.8983 

.9530 

.1332 

.007S 

.900 

20 

30 

.9013 

.9470 

.0938 

.0038 

— 

8 

8 

.8855 

.9650 

.1954 

.0126 

i 

8 

30 

.8983 

.9530 

.1332 

.0078 

20 

30 

wmam 

.9470 

.0938 

.003S 

S 

8 

.9215 

.9650 

.1579 

.0092 

1.0 

8 

30 

.94S6 

.9510 

.0990 

.0049 

.950 

20 

30 

.9498 

.9500 

.0688 

.0022 

8 

8 

.9376 

.9650 

.1578 

.0092 

20.0 

8 

30 

.9486 

.9510 

.0989 

.0049 

20 

30 

.9500 

.06S8 

.0022 

8 

8 

.9846 

.9660 

■0S84 

.0040 

1.0 

S 

30 

.9S9S 

.9510 

.0447 

.0014 

.990 

20 

30 

.9905 

.9450 

.0299 

.0006 

S 

8 

■9S47 

.9660 

-0SS3 

.0040 

20.0 

8 

30 

.9S9S 

.9510 

.0447 

.00  J  4 

20 

30 

.9450 

.0298 

.0006 

S 


A • 

The  above  method  for  deriving  cannot  be  used  to  find  a  lower  confidence  in¬ 
terval  for  P{ X  >  f)  when  we  drop  the  assumption  of  equal  variances.  Consequently, 
we  shall  develop  using  a  different  approach  which  has  greater  potential  for  con¬ 
structing  confidence  intervals  when  variances  are  not  equal.  Let 

K  =  g(X  -  F,  S2)  =  .  (2.8) 


The  Taylor  expansion  of  K  at  {nx  -  /ir,  o2)  ,  using  only  first  order  derivatives,  is  given 
by: 


g(A'  -  ]•',  S*)  =  g[uX  -  fty,  c2)  +  [X  ~  ~  (f*X  ~  Py) 3  -  S 

01  “  /-  , 


3{X  -  Y) 


iPX~PY.e  ) 


2  4-/? 
U<x-Py,c‘)  '  Kt 


(2.9) 


PX~PY  .  X  ~  Y  -  f.u.Y  -  Hr)  /o2  2.  PX~PY  ,  „ 

- - 5 - + - 5 - ( S  -o)  2(_J  +  ‘h 


The  expected  value  and  variance  of  g(A'  -  Y.  S 2)  are  as  follows: 


„  _  Uv —  Uv 

rm  =  f‘s  =  ^ir^ 


Var[X]  = 

=  \  Yarf.Y  -  F)  +  (  )2v'ar(A'2) 

_  1  f  c2  ,_  c2  \  .  (.<';■•  ~  Py)2  2c4 
^2  \  n  ’  in  J  •  ^6  n  +  m  —  2 

_  (  1  ,  l  \  .  {px-py-y 

\  n  +  m  )  "*■  -  2.  .  ,  * 

v  '  2(7  in  -i-  in  —  2 1 


(2.10) 


An  estimator  for  gk,  is  then 


(X  -  Yf 
2S2in  -r  in  -  2) 


K2  . 
2[n  -f  in  -  2)  ’ 


(2.11) 
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I.  Example 

We  illustrate  the  application  of  this  \  '„edure  with  an  example.  We  compute 

A 

the  90%  lower  confidence  limit  of  the  component  reliability,  RLtM  given  :hc  following 
data: 


10 


X:  9.26 

10.  19 

9.  79  ' 

11.  27 

10.06 

9. 22 

9.  75 

8.46  9.  79  9.52 

11.55 

9.  41 

9.  99 

9.  22 

11.22 

9.  89 

9.54 

9.24  8.77  9.86 

10.  03 

10.33 

10.95 

8.33 

10.85 

Y:  8.07 

8.53 

7.7o 

3  4  5 

■f  ’■'> 

6.  36 

6.  82 

8.65  7.00  7.65 

7.  94 

9.  15 

6.  99 

7.  30 

j.  88 

9.56 

9.58 

7.85  S;  91  7.52 

9.  97 

8.  76 

7.  01 

9.96 

7  4.** 

9.  05 

9.54 

10.72  7.86  7.87 

n  =  25,  3F  =  9.860 


m  =  30,  y  =  8.308 


s  = 


V 


Z(x,-x)-~Z(y-) 
'  i=l 


__ 


it  +  m  —  2 


=  0.97S 


K 


X-J 

s 


1.587 


A 

^,0.9  =  ^- 


2(n  +  m  -  2) 


^0.9,  n~m~2 


0.798 


We  note  that  the  point  estimator  of  R, 


A 


c\/2  /  V  V2  ' 


=  <I>(1.122)  =  0.868, 


and  the  computed  90%  lower  confidence  limit  of  R,  RLi0_ ,  =  0.79S  <  0.S6S  as  it  should 
be.  We  cannot  draw  any  conclusion  about  the  accuracy  of  any  confidence  interval 
procedure  from  one  example.  Wo  need  computer  simulations  to  do  '.nib.  The  next  sec¬ 
tion  considers  the  accuracy  of  this  procedure. 
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B.  ACCURACY  OF  THE  PROCEDURE 

1.  'Measures  of  accuracy  and  the  concept  of  computer  simulation 

The  accuracy  of  the  interval  procedure  in  Equation  (2.12)  was  evaluated  in 
terms  of  the  following  four  chaiaueristics: 

•  The  actual  confidence  level  of  the  interval,  i.e.  the  portion  of  times  an  estimated 
limit  will  cover  the  true  reliability. 

•  The  mean  error  between  th»  estimated  limit  and  the  true  reliability,  which  is  de¬ 
noted  as  'mean  error  from  R  in  the  simulation  result  tables. 

The  variance  of  thi  error  between  the  estimated  limit  a  .d  true  reliability,  which  is 
denoted  as  'variance  of  error'  in  the  tables. 

•  The  100(1  —  a),h  percentile  point  of  the  distribution  of 

The  actual  confidence  level  of  an  approximate  confidence  interval  can  be  ac¬ 
cessed  during  a  computer  simulation  which  we  will  discuss  later. 

To  compute  the  100(1  -a),h  percentile  point  of  the  RLA.„  one  only  need  to  ex- 

A 

amine  the  definition  of  Ru  as  a  lower  100(1  -  oc)  %  confidence  limit  for  R.  i.e. 

i  «  <  R)  i-«.  (2.14) 

This  equation  says  that,  R  is  the  100(1  --  a),h  percentile  point  of  the  probability  distrib- 
ution  of  RLi  Thus,  if  we  construct  the  distribution  of  RLA.,  b\  computer  simulation, 
we  should  find  that  the  100(1  —  ct):h  percentile  point  of  our  constructed  distribution  is 

A 

R,  provided  RM_,  is  a  true  100(1  —  cc)%  lower  confidence  limit  for  R.  Let  R'u denote 
the  1 00(1  —  «)!h  percentile  poin;  of  RLA.t-  Then  the  quantity  |  —  /?  |  is  a  measure 

of  the  accuracy  of  the  procedure. 

A 

We  can  construct  the  distribution  of  RLA.,  by  generating  a  large  number,  say 

A 

1000,  of  random  observations  on  R,,,  for  a  given  set  of  parameter  values 
Px.  Pr.  <72i  and  R.  The  100(1  -  a)  empirical  percentile  point  of  the  distribution  of 
Rl.  i-,  is  the  1000(1  -c<)'h  ordered  statistic  of  Rtil_„ ,  R, ooon-.j.  «i-»»  ■ 


2.  Computer  simulation 


I 


a.  Simulation  procedure 

Eighteen  sets  of  values  of  //,  in,  px,  pY,  and  o2  were  chosen  to  perform  the 
simulations.  They  were  selected  in  a  manner  so  that  R-<1>(  )  ,  for  R  =  0.90, 

\  2yj<7  / 

0.95,  0.99.  Thus  when  random  samples  of  X  and  Y  are  generated,  the  reliability 
R  =  I\X  >  l’)  -will.  be  at  designed  vj’ic.  Parameters  were  also  chosen  according  to  the 
following  rules  in  order  to  cover  practical  conditions: 

•  PX>Py 

•  Reliability,  R  :  0.90,  0.95,  0.99 

•  Standard  Deviations,  c:  1.0,  20.0 

•  Sample  sizes,  (n,  m):  'S,  8);  (8,  30):  (20,  30) 


Actual  sets  of  parameters  ar  ■  tabulated  on  Appendix  B. 

Each  set  of  parameters  describes  a  case.  For  each  case,  the  set  of  parameters 
were  used  to  generate  samples  of  size  n  and  m  for  normal  variates  X  and  Y  respectively. 
The  developed  methods  were  used  to  estimate  the  confidence  limit  for  confidence  levels 
of  a=  0.05,  0.1,  and  0.20.  For  each  individual  case,  this  procedure  is  replicated  1000 
times  producing  1000  random  observations  of  The  1000(1  —  a)"1  ordered  statistic 

r 

of  RLj  is  compared  with  the  true  reliability.  The  actual  confidence  levels  is  computed 
by  counting  the  number  of  the  1000  Ru statistics  that  fall  below  the  true  reliability. 
We  also  computed  the  sample  mean  and  variance  of  the  estimation  error  using  the  1000 
generated  values  of  RL^. 

The  simulation  procedure  can  be  summarized  as  follows: 

1)  Generate  n  random  norma!  variates  X,  T'~  N(pxg2)]  m  random  normal  variates 

Y.NAM 

2)  Compute  X,  Y,  S. 

A 

3)  Compute  K, 


4)  Compute 


R 


L,  l-j 


C) 


L.  I-j 


v2 


for  a 


0.05,  0.10,  0.20. 
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5)  Repeat  steps  1  through  4,  for  1000  times. 

AAA  A 

^<2).£(l-»)  »  ^1000(l-aU(l-»)  l  •••!  ^(1000). 


Then  order  the  to  get 

from  the  smallest  to  the  largest. 


A 

6)  Print  ^ioooo-»).i(i-»r 


7)  Print  max  <  #}  ,the  actual  confidence  level  of  the  approximation 

will  be 


1000 


8)  Compute  and  print  the  sample  mean  and  variance  of  the  estimation  error,  to 
measure  the  precision  and  stability  of  the  approximation. 


b.  Simulation  language 

The  cimulations  were  conducted  on  the  N.P.S.  mainframe  IBM  3033  com¬ 
puter.  The  programming  language  used  for  the  simulations  is  VS  FORTRAN  2.  The 
random  number  generator  LLRANII  was  used  to  generate  normal  variates.  The  IMSL 
statistics  function  TIN  was  used  to  compute  the  percentile  of  the  student  t  distribution 
and  ANORDF  was  used  to  compute  the  probability  of  standard  normal  distribution. 
The  FORTRAN  source  code  for  the  simulation  on  the  approximate  procedure  are  at¬ 
tached  in  Appendix  A.  The  FORTRAN  code  for  the  simultaneous  comparison  simu¬ 
lation  of  the  approximate  procedure  and  the  nonparametric  procedure  arc  attached  in 
Appendix  C. 

3.  Analysis  of  simulation  results 

The  simulation  results  are  tabulated  in  Tables  3  through  5.  The  results  shows 
that  is  very  close  to  the  true  reliability  for  every  case  of  simulation,  so  that 

we  have  developed  a  \er\  exact  procedure.  The  procedure  is  more  nearly  exact  for  large 
sample  sizes. 
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Table  3.  ANALYSIS  OF  80%  CONFIDENCE  LIMIT  APPROXIMATION  OF 


R  =  P(X  >  Y)  FOR  EQUAL  VARIANCES  CASE  WITH  EQN.  2.12 


R 

B 

n 

m 

A 

^IOOOO-j).  i( 1-J) 

True  confi¬ 
dence  level 

Mean  error 
from  R 

Variance  of 
error 

.900 

1.0 

8 

8 

.8989 

.8010 

.0840 

.0089 

m 

30 

.9003 

.79S0 

.0628 

.0056 

20 

30 

.9019 

.7880 

.0449 

.0028 

20.0 

8 

8 

.8989 

.8010 

.0840 

.00S9 

8 

30 

.9003 

.7980 

.062S 

.0056 

20 

30 

.9019 

.7880 

.0449 

.0028 

.950 

1.0 

S 

8 

.9498 

.8010 

.0637 

.0052 

S 

30 

.9511 

.7910 

.0450 

.0030 

20 

30 

.9516 

.7850 

.0321 

.0015 

20.0 

8 

8 

.9499 

.SOlO 

.0636 

.0052 

8 

30 

..9511 

.7910 

.0450 

.0030 

20 

30 

.9517 

.7840 

.0320 

.0015 

.990 

1.0 

S 

8 

.9901 

.7970 

.0302 

.0014 

8 

30 

.9906 

.7830 

.0185 

.0006 

20 

30 

.9904 

.7890 

.0130 

.0003 

20.0 

8  'I 

S 

.9901 

.7970 

.0302 

.0014 

8 

30 

.9906 

.7830 

.0185 

.0006 

20 

.7890 

.0130 

.0003 

Table  4;  ANALYSIS  OF  90%  CONFIDENCE  LIMIT  APPROXIMATION  OF 


M00«I-j).  Ul-2 ) 


.8944 


.9050 


.9034 


-S944 


.9050 


.9034 


.9467 


.9529 


.9517 


.9467 


.9530 


.95 


.9SS6 


.9909 


.9904 


.9SS7 


.9909 


.9904 


L  VARIANCES  CASE  WITH  EQN.  2.12 

True  confi¬ 
dence  level 

Mean  error 
from  R 

Variance  of 
error 

.9060 

.1390 

IHBDH 

.SSSO 

.0989 

.006S 

.8880 

.0703 

.0033 

.9060 

.1390 

.0111 

.8880 

.09  S  9 

.0068 

.SSSO 

.0703 

.0033 

.9090 

.1081 

.0074 

.8860 

.0720 

■mu | 

.8890 

.0507 

.0019 

.9090 

.1080 

.0074 

.SS60 

.0719 

.0040 

.SS90 

.0507 

.0019 

.9110 

.0549 

.0025 

.SS10 

.0308 

.0009 

.8950 

.0211 

.0004 

.9110 

.0549 

.0025 

.SSIO 

.0308 

.0009 

.8950 

.0211 

.0004 

Table  5.  ANALYSIS  OF  95%  CONFIDENCE  LIMIT  APPROXIMATION  OF 
R  =  P(X>Y)  FOR  EQUAL  VARIANCES  CASE  WITH  EQN.  2.12 


S  I  8 


8  30 


20  30 


A 

^IOOOtf-i).  /.(1-j) 

True  confi¬ 
dence  level 

Mean  error 
from  R 

Variance  of 
error 

.8884 

.9600 

.1942 

.0128 

.S9S7 

.9520 

.1330 

.0079 

.9015 

.9470 

.0937 

.0038 

.8SS4 

.9600 

.1942 

.0128 

.8987 

.9520 

.1330 

.0079 

.9015 

.9470 

.0937 

.0038 

.9406 

.9620 

.1559 

.0093 

.9489 

.9510 

.0987 

.0049 

.9500 

.9500 

.0686 

.0022 

.9406 

.9620 

.I55S 

.0093 

.9490 

.9510 

.09S6 

.0049 

.9500 

.9490 

.0686 

.0022 

.9865 

.9620 

.0857 

.0040 

.9900 

.9500 

.0444 

■■■ 

.9906 

.9450 

.0297 

.0006 

.9865 

.9620 

.0S57 

.0040 

j  .9900 

.9500 

.0444 

.0014 

.9906 

.9450 

.0297 

t  .0006 

III.  NONPARAMETRIC  DISTRIBUTION  FREE  LOWER  CONFIDENCE 

BOUND  PROCEDURE 

The  distribution  free  confidence  bound  procedure  suggested  by  Govindarajulu  [Ref. 
4:  pp.  229-238],  may  be  applied  to  our  problem.  However  we  should  note  that  the 
Govindarajulu  procedure  requires  large  sample  sizes  and  the  true  reliability  values  re¬ 
moved  from  0  and  1.  For  comparative  purposes  only,  we  evaluated  his  procedure  using 
the  same  computer  simulation  and  analysis  methods  that  were  performed  on  our  dev  el¬ 
oped  procedure.  The  same  data  was  used  to  evaluate  both  procedures  with  emphases 

A 

on  Riooo(i -#).  and  'mean  error'.  The  results  are  displayed  in  Tables  6  and  7. 
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Table  6.  COMPAF  ON  OF  DIFFERENT  PROCEDURES  ON  90%  LOWER 
CONFIDE  IE  LIMIT  ESTIMATION  OF  R  =  P(X>Y),  EQUAL 
VARIANCES  CASE. 


R 

■■ 

Approximate  Estimation 

Nonparametric  Bound 

n 

n 

m 

A 

Mean  Error 

m 

30 

.8998 

.0650 

.8225 

.1289 

1.0 

50 

70 

.8969 

.0419 

.8417 

.0915 

.900 

90 

90 

.9000 

.0322 

.8602 

.0676 

25 

30 

.8998 

.0650 

.8225 

.1289 

20.0 

50 

70 

.8970 

.0418 

.8417 

.0915 

90 

90 

.9000 

.0321 

.8602 

.0676 

25 

30 

.9488 

.0466 

.8545 

.1289 

1.0 

50 

70 

.9484 

.0290 

.8805 

.0913 

.950 

90 

90 

.9499 

.0220 

.8997 

.0677 

25 

30 

.9488 

.0465 

.8545 

.1288 

20.0 

50 

70 

.9485 

.0288 

.8805 

.0912 

.0219 

.8999 

.0677 

25 

30 

.9894 

.0190 

.8718 

.1286 

1.0 

m 

70 

BHHI 

.0107 

.9068 

.0908 

.990 

90 

90 

.9900 

.0079 

.9284 

.0676 

25 

30 

.9895 

.0190 

.8718 

.1286 

20.0 

50 

70 

.0107 

.9068 

.0908 

90 

90 

.9900 

.0079 

.9284 

.0676 
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Table  7.  COMPARISON  OF  DIFFERENT  PROCEDURES  ON  95%  LOWER 
CONFIDENCE  LIMIT  ESTIMATION  OF  R  =  P(X>Y),  EQUAL 
VARIANCES  CASE. 


R 

a 

n 

m 

Approximate  Estimation 

Nonparametric  Bound 

Mean  Error 

A 

■^iooou-j), 

Mean  Error 

.900 

1.0 

25 

30 

,8989 

.0S66 

.7955 

.1653 

50 

70 

.8976 

.0549 

.8245 

.1172 

WEI 

90 

.8994 

.0422 

.8466 

.0S68 

20.0 

25 

30 

.8989 

.0866 

.7955 

.1653 

50 

70 

.8977 

.054S 

.8245 

.1172 

90 

90 

.8995 

.0422 

.8466 

.0868 

.950 

1.0 

25 

30 

.9494 

.0630 

.8235 

.1652 

50 

70 

.9477 

.0384 

.8597 

.1170 

90 

.9490 

.0291 

.8850 

.0S69 

20.0 

25 

30 

.9495 

.0629 

.8235 

.1652 

50 

70 

•947S 

.8597 

.1169 

90 

90 

.9491 

.0290 

.8850 

.0S68 

.990 

1.0 

25 

30 

•9901 

.0267 

.S355 

.1650 

50 

70 

.9895 

.0146 

.8820 

.1165 

90 

90 

.9897 

.0107 

.9102 

.086S 

20.0 

25 

30 

.9901 

.0266 

.8355 

.1649 

■3 

70 

.9895 

■B 

.8S20 

.1165 

90 

90 

.9897 

mmm mm 

.9102 

.0867 
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IV.  APPROXIMATE  INTERVAL  ESTIMATION  PROCEDURE  FOR 
RELIABILITY  R  =  P(X  >  Y) -  UNEQUAL  VARIANCES  CASE 

A.  LOWER  CONFIDENCE  LIMIT  PROCEDURE 

Let  X  denote  component  strength  where  X  ~  X(nxPx)  Let  Y  denote  stress  applied 
to  the  component  where  Y  ~  N(nYo2Y) .  Then 

X  —  Y~  N(nx  —  fJ-y,  o\  +  c2y). 

The  component  reliability  is  defined  as  follows: 

R  =  P[X  >  i1] 


=  p 

X-Y-fa-py)  . 

Rx-Ry 

yj C% +  Cy 

j  2  2 

\l<?x  +  cy  J 

(4.1) 

=  (]>/ 

HX-t‘Y  , 

i  2  1 

v  +  °y 

where  <I>  is  the  standard  normal  cumulative  distribution  function.  Let 


where  and  S-y  denotes  sample  variances  and  A'  and  Y  are  the  respective  sample 
means. 
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For  an  observed  value  of-.a-..-/  the  lower  100(1  —  a)  %  confidence  limit  dL,  for 

V4  +  4 

5,  using  the  general  confidence  interval  method,  is  the  value  of  5  in  Equation  (4.2)  such 
that 


!-«  =  /> 


X  -  Y  ..  x-y 

\l  Sy  +  Sy  yj  fy  +  Sy  . 


(4.4) 


The  probability  statement  in  (4.4)  cannot  be  reduced  to  an  equhalent  statement 

about  a  random  variable  whose  distribution  has  been  tabulated.  This  problem  is  similar 

to  the  Behrens-Fisher  problem  of  finding  a  confidence  interval  for  iix  —  nY  when  both 

variances  a\  and  a\  are  unknown  and  unequal.  B.  L.  Welch  [Ref.  10:  pp.  28-35]  has 

proposed  .an  approximate  confidence  interval  for  /jlx  — ny  in  this  case  using  the  statistic 

—=====- .  Welch  approximated  the  distribution  of  this  statistic  with  a  Student's  t 
Jsiln  +  Slim 

distribution  with  v  degrees  of  freedom.  The  degrees  of  freedom  chosen  by  Welch  is  given 
by 

(Sj-ln  +  S'ylmf 

Sjtln*  Sfrn2 
n—  1  +  »i  -  1 

The  Welch  statistic. is  not  useful  in  our  problem,  because  we  need  to  use  a  statistic  that 

is  a  consistent  estimator  for  ■  ■ . 

v«.f  +  <4-  _  _ 

In  particular  we  choose  to  use  the  statistic  g  =  ——==■ .  The  mean  and  variance 

■Jsi+si 

of  this  statistic  are  approximated  and  the  distribution  of  is  approximated  with 

a  Student's  t  distribution.  The  desired  confidence  interval  is  constructed  using  the  ap¬ 
proximated  Student's  t  distribution.  The  analysis  in  the  following  paragraph  serves  onl> 
as  a  means  to  find  a  plausible  expression  for  the  degrees  of  freedom  of  the  approximated 
Student's  distribution. 
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(x  —y) 

We  begin  as  in  equation  (4.4),  where  — p.=_.  .  denotes  a  value  constructed  from 

V4  +  4 

observed  data.  The  upper  case  version  of  the  same  expression  denotes  a  random  vari¬ 
able. 


1  -«  =  P 


x  -  y 


X  -y 


L  -\/ Sy  +  Sy  V sy  +  *y  J 

A’  -  y  -  (l‘x  -  My)-+  (.«,Y  ~  /'?■) 


=  P 


=  P\ 


\[  7«  +  Gylm 


X  -y 


\‘ 


j  sj-jn  +  Sylm 


isl  +  Sy 


V  °.Y 


Stf/tt  +  Syjm 


f~2  2 

\>s,yr  SY 


Z  + 


y  c}:ln  -j-  try/'" 
t>X-PY 


■J Gxl II  +  Cyjm 


x -y 


SX  +  Sy 

'l  “f 

<7 yin  -r  a  yin 


sxln  4“  SyjtH 

TT3 

w  y  sy 


-Jsxln  +  sft 


ylm 


(4-5) 


where  o\,  a\,  4,  and  denote  constants.  The  last  probability  in  Equation  (4.5)  suggests 
.  a  noncentral  t  distribution.  Consequently  the  random  variable 

I  [$x I'1  4~  Syl>n)(Sx  ~  Sy) 

V  (pxl'1  +  c  yl'»)(Sy  +  Sy) 


in  the  denominator  should  be  a  random  variable  of  the  form 

Sy/n  +  Sylm  Sy  -i -  Sy  2 

2,  ,  2 1  2.2  v~  Xv‘ 

<7 yin  4-  <7  ypn  Sx  T  Sy 


That  is, 


Since  Var (rf)  =  2v  ,  Var(Sj)  =  .  and  Var(Sj)  = 


2v  =  v 


2/  SX In  +  Sylm  V 

„2  ,  .2 

~  5>' 


2<7y  2(7  y 

{<72yln  +  <72ylm)2  V  ~+  -  1 


Thus, 
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fyjtjr 

4l'l  +  Sylm 

4  +  4 


(<s\l»  +  (r  j »/»i)2 
<4/(«  -  1)  +-<7y/(»l  -  1) 

i4'ln  +  syl»:)2 


sX>n  +  *}’/">  /  4/(«  -  1)  +  iy/(w  “  1) 

_ (4+4? 


sx/(»-V  +  sy/(m-  1) 


(4-6) 


In  summarv,  if  we  fit  a  t  distribution  to  the  distribution  of -g  ,  the  expression  in 
the  Equation  (4.6)  might  be  a  plausible  random  value  for  the  degrees  of  freedom.  The 
results  of  computer  simulationswill  indicate  the  accuracy  of  this  approximation. 


In  order  to  formulate  the  lower  confidence  interval  statement. for  — 


f'r  -  t‘y 


_  +  a? 


,  we  first 


find  the  mean,  E[g]  =  nv  and  variance,  Var[g]  =  c\,  for  g  = 


X  -  Y 


■.  We  then  ap¬ 


proximate  the  distribution  of-— — £  with  a  central  Student's  t  distribution  with  v  degrees 
of  freedom,,  where  v  is  given  in  equation  (4.6).  The  lower  confidence  limit,  frSL  ,_s  for 
gs  will  yield  the  corresponding  lower  confidence  limit  for  reliability 

4 1-,= 

We  now  proceed  to  find  n.  and  c\.  The  Taylor  expansion  of  g  at  (nx  -  nr,  o\  +  o])  ,  us¬ 
ing  only  first  order  derivatives  is  given  by: 


g(X  -  Y,  Sx  +  Sy)  =  s(nx  -  I*yfx  +  a  y)  +  IX  -  Y  -  (ftx  ~  ^y)3  .-^S  - : 

d(\  -  Y) 


°X~cy'I 


Vx  ~  Vy  ,  X  -  Y  -  [fix~  ft y) 


(ftx  ~  My,  °x  ~  °>')  Rf 


(4.7) 


■JGX  + 


<Py 


-  -  +  Sy  -  [<?x-+  Cy)l  — r - fTF  +  Ri 

yJox  +  <fr  *(eZ  +  <ryr 


where  R,  includes  terms  that  converge  to  0  at  the  same  rate  as  ,  as  n  and 

m  become  large. 

The  expected  value  and  variance  of  g  are 
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(4.8), 


1.  Example 

Let  x,  denote  the  strength  of  solid  missile  motor  chambers  that  were  pressurized 
until  they  burst,  /  =  1,  2,  25.  Let»  denote  the  maximum  pressure  observed  on  30 

solid  missile  motors  that  use  this  type  of  chamber.  X  and  Y  are  assumed  to  have  normal 
distributions  with  unequal  variances.  The  coded  X  and  Y  data  arc: 
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292. 65 

301.  86 

297.  86 

312.  67 

300. 

61 

292. 18 

297.55 

284. 64 

297.  93 

295. 22 

315.  49 

294. 

10 

299. 85 

292.  24 

312.25 

298. 87 

295.41 

292.  41 

287. 

71 

298.57 

300. 34 

303. 28 

309.52 

283.  32 

308. 52 

242.52 

261.00 

229.  37 

26l.:82 

248. 

36 

174. 04 

192.50 

265. 64 

199. 70 

225.52 

237-43 

285. 

66 

199. 20 

203.  65 

274. 93 

302.  13 

302.  85 

233. 59 

275. 

95 

220.46 

318.  34 

269.86 

200. 10 

318.  19 

216. 69 

281. 

78 

301.  20 

348. 56 

234. 15 

234. 45 

The  corresponding  90%  lower  confidence  limit  for  component  reliability, 
R  =  P( X  >  Y)  is  computed  as  follows: 

n  =  25,  x  =  298.C0,  s%  -  70.25 

m  =  30,  jT  =  251.99,  sy=  1873.79 


v= - - =31 

4i(n  - 1) + 4i(,n  - 1) 


A 


sfxln  +  Syjm 

2  ,  2 
SX  +  Sy 


(*  -  y )2  f  4  .  4  \ 
2(shsfr  V"-1  “-1  ) 


1/2 


0.23 


s=- 


x  -y 


\!SX  +  Sy 


■=1.06 


*L,  0.9  -S~  V<3-9. 31  =  °-76 


fy.,  0.9  =  ®{6L,  O.o)  =  0-7S 
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The  corresponding  90°, o  Govindarajulu  nonparametric  confidence  limit  is  0.69.  This  is 
not  surprising  because  our  procedure  uses  more  information  about  the  distributions  of 
X  and  V.  Intervals  ..stimated  b\  nonparametric  procedures  are  usually  wider  than  those 
estimated  by  parametric  procedures.  The  amount  of  the  difference  is  somewhat  sur¬ 
prising  for  these  sample  sizes  of  25  and  30. 

B.  ACCURACY  OF  THE  PROCEDURE 
1.  Computer  simulation 
a.  Si/nulaticn  procedure. 

Computer  simulations  were  used  to  determine  its  accuracy  for  specific  sets 
values  of  n,  m,  px,  nr,  o\,  and  n\. 

Parameters  are  chosen  in  a  way  to  cover  practical  conditions.  The  sets  of 
parameters  are  as  follows: 

•  Hx>Vr 

•  Reliability,  R  :  0.90,  0.95,  0.99 

•  Standard  Deviations 

ox-  1  0  10.0 
oy\  2:0,  40.0 

•  Sample  sizes,  (n.  m):  (10,  20);  (25,  35):  (75,  50) 

•  Confidence  leva,  :  0.05,  0.10,  0.20  . 

Actual  sets  of  parameters  are  tabulated  in  Appendix  E. 

Simulation  procedure  is  summarized  as  follows: 

1)  Generate  n  random  normal  variates  X,  X  ~  X{nxp7x)\  m  random  normal  variates 
Y,  X{pyg\). 

2)  Compute  X,  Y,  S2X.  Sj-. 

* 

3)  Compute  os5 L  }_,. 

4)  Compute  RLi1_ ,  =  cl>(dtil_,),  for  a  =  0.05,  0.10,  0.20. 

* 

5)  Repeat  steps  1  through  4,  for  1000  times.  Then  order  the  RLJ.3  to  get 

kihio-,*  Raua-v  from  the  smallest  to  the  largest. 


6)  Print  R1000(,_s)ii(,_,,. 

7)  Print  max^.^,..,:  #01ii(1_a)  <  R }  ,the  actual  confidence  level  of  theapproximation 
Wlli  be  1000' 

8)  Compute  and  print  the  sample  mean  and  variance  of  the  estimation  error. 

b.  Simulation  language 

The  simulations  were  conducted  on  the  N.P.S.  mainframe  IBM  3033  com¬ 
puter.  The  programming  language  used  for  the  simulations  is  VS  FORTRAN  2.  The 
random  number  generator  LLRAXII  was  used  to  generate  normal  variates.  The  IMSL 
statistics  function  TIN  was  used  to  compute  the  percentile  of  Student's  t  distribution, 
and  ANORDF  was  used  to  compute  the  probability  of  standard  normal  distribution. 
The  FORTRAN  source  code  for  the  simulation  of  the  developed  approximate  procedure 
are  attached  in  Appendix  D.  The  FORTRAN  code  for  simultaneous  comparison  be¬ 
tween  the  approximate  procedure  and  nohparametric  procedure  are  attached  in  Appen¬ 
dix  F. 

2.  Analysis  of  simulation  results 

The  simulation  results  are  tabulated  in  Tables  S  through  10.  The  results  show 
the  developed  approximate  procedure  is  quite  accurate,  because  for  every  case  the 

A 

100(1  —  «)<h  percentile  point  of  are  ah  very  close  to  the  true  reliability.  Also,  the 
'mean  error  from  R'  and  'Variance  of  error'  reduce  rapidly  with  the  increases  of  sample 
size.  Comparison  simulations  were  run  for  the  approximate  procedure  versus  the 
Govindarajulu  nonparametric  procedure.  Results  are  tabulated  in  table  II  and  12. 
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Table  8.  ANALYSIS  OF  80%  CONFIDENCE  LIMIT  APPROXIMATION  OF 


R  =  P(X>  Y)  FOR  UNEQUAL  VARIANCES  CASE 


R 

n 

n 

m 

A 

■^10000-1),  ;.(t-a) 

True  confi¬ 
dence  level 

Mean  error 
from  R 

Variance  of 
error 

10 

20 

.900S 

.7960 

.0552 

.0044 

1.0 

2.0 

25 

35 

.8996 

.8030 

.0388 

.0021 

:900 

75 

50 

.S997 

.S010 

.029$ 

.0011 

10 

20 

.9019 

.7910 

.0530 

.0042 

10.0 

40.0 

25 

35 

.9012 

.7900 

.0392 

.0022 

75 

50 

.9000 

.SOOO 

.0317 

.0013 

10 

20 

.9500 

.S000 

.0409 

.0024 

1.0 

2.0 

25 

35 

.9502 

.79S0 

.0278 

.0011 

.950 

m 

.9496 

.8080 

.0211 

.0006 

.9513 

.7840 

.0401 

.0024 

10.0 

40.0 

OS 

35 

.9514 

.7S20 

.02S7 

■  .0012 

75 

50 

.94S7 

m 

.0229 

.0007 

10 

20 

.9902 

.7930 

.0180 

.0005 

1.0 

2.0 

25 

35 

.9903 

.7890 

.0112 

.0002 

.990 

75 

50 

.9S9S 

.8110 

.0082 

.0001 

10 

.9906 

.0186 

.0006 

10.0 

40.0 

25 

35 

.9907 

.7750 

.0122 

.0002 

m 

50 

.9899 

.S030 

.0093 

.0001 

Table  9.  ANALYSIS  OF  90%  CONFIDENCE  LIMIT  APPROXIMATION  OF 


R  =  P(X>  Y)  FOR  UNEQUAL  VARIANCES  CASE 


R 

n 

°Y 

n 

m 

A 

•^IOOW-j).  i(l-s) 

True  confi¬ 
dence  level 

Mean  error 
from  R 

Variance  of 
error 

10 

20 

-900S 

.8970 

.0893 

.0052 

1.0 

2.0 

25 

35 

.8993 

.9010 

.0611 

.0024 

.900 

75 

50 

.8990 

.9070 

.0465 

.0013 

10 

20 

.9022 

.8890 

.0869 

.0050 

10.0 

40.0 

m 

35 

.9005 

.8970 

.0624 

.0025 

75 

50 

.8995 

.9020 

.0499 

.0015 

10 

20 

.9497 

.9010 

.0671 

.0032 

1.0 

2.0 

25 

35 

.9502 

.8990 

.0014 

.950 

El 

50 

.9493 

.9030 

.0332 

.0007 

II 

10 

20 

.9510 

.8960 

.0667 

.0032 

■top! 

40.0 

25 

35 

.9494 

.9020 

.0462 

.0015 

■ 

75 

50 

.9494 

.9040 

.0364 

■Eiiai 

10 

20 

.9899 

.9060 

.030S 

.0009 

1.0 

2.0 

25 

35 

.9S99 

.9010 

.0184 

.0003 

.990 

75 

50 

.9S99 

.9050 

.0132 

.0001 

10 

20 

.9901 

.S9S0 

.0321 

.0009 

10.0 

40.0 

25 

35 

.9904 

.S930 

.0203 

SKKEBS 

75 

50 

.9898 

.9060 

.0152 

.0002 

Table  10.  ANALYSIS  OF  95%  CONFIDENCE  LIMIT  APPROXIMATION  OF 


R  =  P(X>  Y)  FOR  UNEQUAL  VARIANCES  CASE 


R 

m 

n 

in 

A 

^lOOOil-i).  1(1-1) 

True  confi¬ 
dence  level 

Mean  error 
from  R 

Variance  of 
error 

10 

20 

.8958 

.9560 

.1218 

.0059 

1.0 

2.0 

25 

35 

-S9SS 

.9530 

.0816 

.0027 

.900 

75 

.S99S 

.9500 

.0615 

.0014 

10 

.9470 

.1196 

.0057 

10.0 

40.0 

25 

BM 

-S9S0 

.9540 

.0838 

.002S 

75 

50 

.9520 

.0664 

.0017 

10 

20 

.9461 

.9550 

.0934 

.0039 

1.0 

2.0 

E9 

ES 

.94S2 

.9590 

.0600 

.0016 

.950 

75 

50 

.9500 

.9480 

.0445 

.0008 

10 

20 

.9517 

.9450 

.0938 

.0039 

10.0 

ii. 

o 

o 

25 

35 

.9470 

.9550 

.0632 

.0018 

75 

50 

.9500 

.9500 

.0492 

.0010 

10 

20 

.9SSS 

.9550 

.0452 

.0013 

1.0 

2.0 

25 

35 

.9S9 1 

.9610 

.0259 

.0004 

.990 

75 

50 

.9S99 

.9520 

.0IS2 

1  .0002 

10 

2ii 

.9903 

.9420 

.0478 

.0014 

10.0 

40.0 

25 

35 

.9892 

.9550 

.02SS 

■■ 

75 

|  .9490 

.0212 

.0002 
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Table  11.  COMPARISON  OF  DIFFERENT  PROCEDURES  ON  90%  LOWER 
CONFIDENCE  LIMIT  ESTIMATION  OF  R  =  P(X  >  Y),  UNEQUAL 
VARIANCES  CASE 


R 

HI 

HI 

II 

Approximate  Estimation 

Nonparametric  Bound 

m 

H 

n 

i 

A 

-^kxxxi-i).  u  i-j) 

Mean  error 

A 

•^ioocki-j). 

Mean  error 

10 

15 

.9013 

.1004 

.7707 

.2018 

1.0 

2.0 

70 

35 

.9003 

.0564 

.8370 

.1098 

.900 

90 

90 

.9009 

.0333 

.8623 

.0673 

10 

15 

.9056 

.1010 

.7840 

.2017 

10.0 

o 

o 

70 

35 

.9010 

.0613 

.8452 

.1094 

90 

90 

.9001 

.0349 

.8684 

.0672 

10 

15 

.9523 

.0767 

.7974 

.2015 

1.0 

2.0 

70 

35 

.9506 ; 

.0410 

.8729 

.1091 

.950 

90 

90 

.9507 

.0233 

.9027 

.0675 

10 

15 

.9550 

.0790 

.7974 

.2014 

10.0 

40.0 

70 

35 

■Em 

.0456 

.SS07 

.1090 

90 

90 

— 

.0250 

.9067 

.0672 

10 

15 

.9915 

.0370 

.7974 

.2021 

1.0 

2.0 

70 

35 

.9906 

.0172 

.8913 

.1088 

.990 

90 

90 

.9901 

.ooss 

.9301 

.0676 

10 

15 

.9921 

.0402 

.7974 

.2021 

10.0 

K 

70 

35 

.9904 

.0200 

.8917 

.1087 

1 

90 

90 

11 

.0098 

.9320 

.0675 
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Table  12.  COMPARISON  OF  DIFFERENT  PROCEDURES  ON  95%  LOWER 
CONFIDENCE  LIMIT  ESTIMATION  OF  R  =  P(X>  Y),  UNEQUAL 
VARIANCES  CASE 


R 

ax 

m 

n 

m 

Approximate  Estimation 

Nonparametric  Bound 

A 

Mean  Error 

A 

R|0«Kl-3>.  //!—»> 

Mean  Error 

.900 

1.0 

2.0 

10 

15 

.8988 

.1388 

.7266 

.2593 

70 

35 

.8977 

.0750 

.8194 

.1405 

90 

90 

.9014 

.0439 

.8529 

.0865 

10.0 

40.0 

10 

15; 

.8975 

.1414 

.7399 

.2592 

70 

35 

.8950 

.0823 

.8275 

.1401 

90 

90 

.9005 

.0463 

.8575 

.0S63 

.950 

1.0 

2.0 

10 

15 

.9506 

.1086 

.7399 

.2590 

70 

35 

.9477 

.0553 

.8483 

.1398 

90 

90 

.9503 

.0310 

.SS90 

.0866 

10.0 

40.0 

10 

15 

.9503 

■H 

.7399 

.2589 

70 

35 

.9468 

.0622 

.8557 

.1397 

90 

90 

.9509 

.0334 

.8927 

.0864 

.990 

1.0 

2.0 

10 

li 

.9S96 

.0557 

.7399 

.2595 

70 

ES 

.9898 

.0240 

.S610 

.1395 

.9903 

.0119 

.9117 

.0S67 

10.0 

40.0 

10 

15 

.9907 

.0617 

.7399 

.2595 

70 

35 

.9893 

.0285 

.S610 

.1395 

90 

90 

.9901 

.0134 

.9132 

.0S67 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  lower  confidence  limit  estimation  procedures  developed  in  this  thesis  for  equal 
variances  case  as  well  as  for  unequal  variances  cases,  are  very  accurate.  These  proce¬ 
dures  are  simple  to  evaluate  and  require  only  the  use  of  central  Student's  t  tables,  in 
contrast  to  the  existing  parametric  procedures  of  this  type  which  require  extensive  use 
of  noncentral  Student's  t  tables. 

Although  these  procedures  are  developed  for  lower  confidence  limits,  upper  or  two- 
sided  confidence  limits  for  the  reliability  are  readily  obtained. 
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APPENDIX  A.  FORTRAN  CODE  FOR  INTERVAL  ESTIMATION 

PROCEDURE  -  NORMAL  EQUAL  VARIA  NCES  CASE 

PROGRAM  EQSIGM 

*  Vr  it  itit  itit  it  it  it  it  Vo  c  j  c  it  it  it  it  ititit  *  iV  5*c  Vf  Vr  it  Vc  VcVnV  VcsV  ^Wc  itit  7WrVr  Vr  Vc  Vr*  VoVVr  * j'c  }Wc  Vc  ititit  ititititititit  VoV  Vc  iWr  itit 

*  * 

*  THIS  PROGRAM  IS  TO  VALIDATE  THE  LOWER  CONFIDENCE  BOUND  APPROXI-  * 

*  MATION  PROCEDURE  FOR  P(  X  >  Y  ),  WHERE  X,  Y  ARE  NORMALLY  DIS-  * 

*  TRIBUTED  WITH  UNKNOWN  MEANS  AND  A  UNKNOWN  BUT  EQUAL  VARIANCE.  * 

*  * 

*  VARIABLES  DESCRIPTION: 


it 

ALPHA 

- 

NOMINAL  CONFIDENCE  LEVEL 

if 

* 

ANORDF 

- 

IMSL  FUNCTION  FOR  NORMAL  PROBABILITY 

if 

it 

CASE 

- 

NUMBER  OF  TEST  PARAMETER  SETS 

if 

ER 

- 

ERROR  BETWEEN  LIMIT  AND  TRUE  RELIABILITY 

if 

* 

ERBAR 

- 

AVERAGE  OF  ER 

if 

ERSSQ 

- 

SUM  OF  SQUARES  OF  ER 

it 

V? 

ERSUM 

- 

SUM  OF  ER 

it 

* 

ERSV 

- 

SAMPLE  VARIANCE  OF  ER 

i: 

* 

CLOSE 

- 

INDEX  OF  THE  CLOSEST  ESTIMATE 

it 

* 

DELTA 

- 

NONCENTRALITY  OF  T  DISTRIBUTION 

it 

if 

DF 

- 

DEGREE  OF  FREEDOM  OF  T  DISTRIBUTION 

it 

* 

K 

- 

STATISTIC  TO  ESTIMATE  DELTA 

it 

* 

LNORM 

- 

RANDOM  NUMBER  GENERATOR  FOR  NORMAL  VARIATES 

it 

* 

M 

- 

SAMPLE  SIZE  OF  Y  RANDOM  VARIABLE 

it 

it 

MUX 

- 

POPULATION  MEAN  OF  X  RANDOM  VARIABLE 

it 

it 

MUY 

- 

POPULATION  MEAN  OF  Y  RANDOM  VARIABLE 

it 

it 

N 

- 

SAMPLE  SIZE  OF  X  RANDOM  VARIABLE 

it 

* 

R 

- 

REAL  RELIABILITY 

it 

* 

REP 

- 

REPETITION  OF  SIMULATION 

it 

tfc 

RLHAT 

- 

LOWER  CONFIDENCE  LIMIT  OF  RELIABILITY 

it 

if 

SIGMAX 

- 

POPULATION  STANDARD  DEVIATION  OF  X 

it 

if 

SIGMAY 

- 

POPULATION  STANDARD  DEVIATION  OF  Y 

it 

if 

SP 

- 

POOLED  SAMPLE  VARIANCE  OF  X  AND  Y 

it 

if 

SUMSQX 

- 

SUM  OF  SQUARES  OF  X 

it 

if 

SUMSQY 

- 

SUM  OF  SQUARES  OF  Y 

it 

if 

SUMX 

- 

SUM  OF  X 

it 

if 

SUMY 

- 

SUM  OF  Y 

it 

if 

TIN 

- 

IMSL  FUNCTION  TO  COMPUTE  PERCENTILE  OF  T  DIST. 

it 

if 

XBAR 

- 

AVERAGE  OF  X 

it 

if 

YBAR 

- 

AVERAGE  OF  Y 

it 

* 


*  REQUIRED  EXTERNAL  FUNCTIONS:  LNORM,  TIN,  ANORDF  * 

it  it  :rVc  ititit  it  itit  •jV  ?  Wc  ;>c  Vo'c  it  it  Vr  Vc  itit  Vr  Vc  ■}  f  VrsV  it  it  it  VoWr  ?'c  itit  Vc  ■»  r  Vc  Vc  s'oWrj’c  itititit  s‘nV  ?’c  Vrs  r  Vr  tV  it  Vc  VcVc  Vc  it  :’f  Vr  x  #.  x 


INTEGER  REP,  CASE 
REAL  ALPHA,  TOO 
PARAMETER( ALPH A=0 . 05) 
PARAMETER( TO0=2 .  0) 
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PARAMETER  RE P=1000) 

PARAMETER ( CASE=18 ) 

* 

INTEGER  I,  II,  J,  U,  V,  XSEED,  YSEED,  N,  M,  CLOSE 

REAL  MUX,  MUY,  XBAR,  YBAR,  SIGMAX,  SIGMAY,  DF,  R,  RLHAT(REP) , 

+  K,  DELTA,  XX(100),  YY(IOO),  X,  Y,  SUMSQX,  SUMSQY,  SUMX, 

+  SUMY,  TIN,  ANORDF,  RN,  RM,  DIFF,  SP,  TEMP,  ER,  ERSUM, 

+  ERSSQ,  ERBAR,  ERSV 

* 

CALL  EXCMS( 'FILEDEF  12  DISK  SETUP1  DATA  Al‘) 

CALL  EXCMSC FILEDEF  18  DISK  OPT1  DATA  Al') 

* 

DO  2000  1=1,  CASE 

READ  (12,2200)  XSEED,  YSEED,  MUX,  MUY,  SIGMAX,  SIGMAY,  N,  M,  R 
DF  =  REAL(N+M-2) 

RN  =  REAL(N) 

RM  =  REAL(M) 

ERSUM  =  0. 0 
ERSSQ  =  0.  0 
DO  1000  J=l,  REP 

CALL  LNORMC XSEED,  XX,  N,  2,  0) 

CALL  LNORM( YSEED,  YY,  M,  2,  0) 

SUMSQX  =  0. 0 
SUMSQY  =  0. 0 
SUMX  =  0.  0 
SUMY  =  0.  0 

* 

*  <  TRANSFORM  X,  Y  TO  DESIRED  PROPERTIES  > 

Vr 

DO  200  U=  1,  N 

X  =  XX(U)  *  SIGMAX  +  MUX 
SUMSQX=  SUMSQX  +  X  *  X 
SUMX  =  SUMX  +  X 

200  CONTINUE 

XBAR=  SUMX  /  RN 
DO  300  V=l,  M 

Y  =  YY(V)  *  SIGMAY  +  MUY 
SUMSQY  =  SUMSQY  +  Y  *  Y 
SUMY  =  SUMY  +  Y 

300  CONTINUE 

Vc 

*  <  COMPUTE  CONFIDENCE  LIMIT  FOR  RELIABILITY  > 

YBAR  =  SUMY  /  RM 

SP  =  SQRT(  (SUMSQX  -  RN*XBAR*XBAR  +  SUMSQY  -  RM*YBAR*YBAR) 

+  /  DF) 

K  =  MAX(  (XBAR  -  YBAR)  /  SP,  0.0  ) 

DELTA  =  K  -  SQRT(  ( RN+RM) / ( RN*RM)  +  K*K  /  (2. 0*(RN+RM-2.  0) ) ) 
+  *  TIN(1.  0- ALPHA, DF) 

RLHAT(J)  =  ANORDF ( DE  LTA / S  QRT ( TWO ) ) 

■jV 

*  <  COMPUTE  THE  MEAN  AND  VARIANCE  OF  ER  > 

ER  =  R  -  RLHAT(J) 

ERSUM  =  ERSUM  +  ER 
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ERSSQ  =  ERSSQ  +  ER  *  ER 
1000  CONTINUE 

ERBAR  =  ERSUM  /  REAL(REP) 

ERSV  =  (  ERSSQ  -  REAE(REP)  *  ERBAR  *  ERBAR  )  /  REAL(REP-1) 

* 

*  <  SORT  CONFIDENCE  LIMITS  IN  ASCENDING  ORDER  > 


DIFF  =  2.  0 
DO  1800  H=l,  REP 

DO  1500  J=I1+1,  REP 

IF  (RLHAT(J)  .LT.  RLHAT('Il) )  THEN 
TEMP  =  RLHAT(Il) 

RLHAT(Il)  =  RLHAT(J) 

RLHAT(J)  =  TEMP 
ENDIF 

1500  CONTINUE 
#v 

*  *  <  FIND  THE  CLOSEST  CONFIDENCE  LIMIT  ESTIMATE  > 


1800 

2000 

2100 

2200 


IF  ((CR-RLHAT(Il))  .  GE.  0. IE-6)  .AND. 

+  ((R-RLHAT(Il))  .LE.  DIFF))  THEN 
DIFF  =  R  -  RLHAT(Il) 

CLOSE  =11 
ENDIF 
CONTINUE 

WRITE  (18,2100)  I,  MUX,  N,  SIGMAX,  MUY,  M,  SIGMAY,  R, 

+  RLHAT(NINT(REAL(REP)*(1.  0-ALPHA))),  RLHAT( CLOSE) , 

+  REAL( CLOSE)/ 1000.  0,  ERBAR,  ERSV 

CONTINUE 

FORMAT(’ SIMULATION:  ’,12,/, 'MUX:  1 ,F5. 1,T16, ' N:  ',12, 

+  T35,1 SIGMAX:  ' ,F4. 1,/,'MUY:  1 ,F7.  3,T16, 'M:  ’,12, 

+  T35, 'SIGMAY:  ' ,F4.  1,/ , 'TRUE  R  :  ',F7.5, 

+  T35, 'RLHAT:  ' ,F7.  5 ,/, 'CLOSEST  RLHAT:  ',F7.5, 

+  T35, 'TRUE  CONFIDENCE  LEVEL:  ’,F5.3,/, 

+  ' MEAN  ERROR  WIDTH:  ' , F7 . 5 , / , 

+  'VARIANCE  OF  ERROR:  ’ ,F7. 5 , ///) 

FORMAT  (15, IX, 15, IX, F5.  1,1X,F7.  3,1X,F4.  1,1X,F4.  1, IX, 12 , IX, 12 , IX, 
+F5.  3) 

STOP 

END 
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APPENDIX  B.  SIMULATION  PARAMETER  SETS  FOR  EQUAL 


VARIANCES  CASE 


FILE: 

SETUP 1 

.  DATA 

(FOR  APPROXIMATE  PROCEDURE) 

16807 

93943 

10.  0 

8.  187 

1.  0 

1.  0 

8 

8 

.900 

1 

16807 

93943 

10.  0 

8.  187 

1.  0 

1.  0 

8 

30 

.900 

2 

16807 

93943 

10.  0 

8.  187 

1.0 

1.  0 

20 

30 

.900 

3 

16807 

93943 

100.  0 

63. 740 

20.  0 

20.  0 

8 

8 

.  900 

4 

16807 

93943 

100.  0 

63. 740 

20.  0 

20.  0 

8 

30 

.  900 

5 

16807 

93943 

100.  0 

63. 740 

20.  0 

20.  0 

20 

30 

.  900 

6 

16807 

93943 

10- 0 

7.  674 

1.  0 

1.  0 

8 

8 

.  950 

7 

16807 

93943 

10.  0 

7.  674 

1.0 

1.  0 

8 

30 

.  950 

8 

16807 

93943 

10.  0 

7.  674 

1.  0 

1.  0 

20 

30 

.950 

9 

16807 

93943 

100.  0 

53.472 

20.0 

20.0 

8 

8 

.950 

10 

16807 

93943 

100.  0 

53.472 

20.0 

20.  0 

8 

30 

.950 

11 

16807 

93943 

100.  0 

53.472 

20.  0 

20.  0 

20 

30 

.950 

12 

16807 

93943 

10.  0 

6.  711 

1.  0 

1.  0 

8 

8 

.990 

13 

16807 

93943 

10.  0 

6.  711 

1.0 

1.  0 

8 

30 

.990 

14 

16807 

93943 

10.  0 

6.  711 

1.0 

1.  0 

20 

30 

.990 

15 

16807 

93943 

100.  0 

34. 211 

20.0 

20.  0 

8 

8 

.990 

16 

16807 

93943 

100.  0 

34. 211 

20.  0 

20.0 

8 

30 

.  990 

17 

16807 

93943 

100.  0 

34. 211 

20.  0 

20.  0 

20 

30 

.990 

18 

X 

Y 

X 

Y 

X 

Y 

SEED 

SEED 

MEAN 

MEAN 

STDV 

STDV 

N 

M 

R 

FILE:  SEQ4  DATA  (FOR  COMPARISON  SIMULATIONS) 


16807 

93943 

10.  0 

8.  187 

1.  0 

1.0 

25 

30 

.  900 

1 

16807 

93943 

10.  0 

8.  187 

1,0 

1.  0 

50 

70 

.  900 

2 

16807 

93943 

10.  0. 

8.  187 

1.0 

1.  0 

90 

90 

.  900 

3 

16807 

93943 

100.  0 

63.  740 

20.  0 

20.0 

25 

30 

.  900 

4 

16807 

93943 

100.  0 

63.  740 

20.  0 

20.0 

50 

70 

.  900 

5 

16807 

93943 

100.  0 

63.  740 

20.  0 

20.0 

90 

90 

.  900 

6 

16807 

93943 

10.  0 

7.  674 

1.  0 

1.  0 

25 

30 

.  950 

7 

16807 

93943 

10.  0 

7.  674 

i.  0 

1.  0 

50 

70 

.  950 

8 

16807 

93943 

10.  0 

7.  674 

1.  0 

1.  0 

90 

90 

.  950 

9 

16807 

93943 

100.  0 

53.  472 

20.  0 

20.  0 

25 

30 

.  950 

10 

16807 

93943 

100.  0 

53.472 

20.  0 

20.  0 

50 

70 

.  950 

11 

16807 

93943 

100.  0 

53.  472 

20.  0 

20.0 

90 

90 

.  950 

12 

16807 

93943 

10.  0 

6.  711 

1.  0 

1.  0 

25 

30 

.  990 

13 

16807 

93943 

10.  0 

6.  711 

1.  0 

1.  0 

50 

70 

.  990 

14 

16807 

93943 

10.  0 

6.  711 

1.  0 

1.0 

90 

90 

.  990 

15 

16807 

93943 

100.  0 

34.  211 

20.  0 

20.  0 

25 

30 

.  990 

16 

16807 

93943 

100.  0 

34.  211 

20.  0 

20.  0 

50 

70 

.  990 

17 

16807 

93943 

100.  0 

34. 211 

20.  0 

20.0 

90 

90 

.  990 

18 

15 

15 

F5.  1 

F7.  3 

F4.  1 

F4.  1 

12 

12 

F5.  3 

X 

Y 

X 

Y 

X 

Y 

N 

M 

R 

SEED 

SEDD 

MEAN 

MEAN 

STDV 

STDV 
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APPENDIX  C.  FORTRAN  CODE  FOR  SIMULTANEOUS  COMPARISON 
SIMULATION  OF  APPROXIMATE  PROCEDURE  VS. 
NONPARAMETRIC  PROCEDURE  -  EQUAL  VARIANCES 

PROGRAM  COMEQ 

*  THIS  PROGRAM  IS  TO  COMPUTE  THE  LOWER  CONFIDENCE  LIMIT  OF  * 

*  R  =  P(  X  >  Y  ) ,  WHERE  X,  Y  ARE  INDEPENDENTLY  NORMALLY  DIS-  * 

*  TRIBUTED  WITH  UNKNOWN  MEANS  AND  A  UNKNOWN  BUT  EQUAL  VAR I AN-  * 

*  CES  WITH  APPROXIMATE  PROCEDURE  AND  WITH  NONPARAMETRIC  PRO-  * 

*  CEDURE  * 

*  * 

*  VARIABLES  DESCRIPTION:  * 


ALPHA  -  NOMINAL  CONFIDENCE  LEVEL 

ANORDF  -  IMSL  FUNCTION  FOR  NORMAL  PROBABILITY 

ANORIN  -  IMSL  FUNCTION  FOR  INVERSE  NORMAL  CDF 

CASE  -  NUMBER  OF  TEST  PARAMETER  SETS 

BIGU  -  MANN-WHITNEY  STATISTIC 

CLOSE  -  INDEX  OF  THE  CLOSEST  ESTIMAT  (APPROXI  PROCEDURE) 

CLOSEN  -  INDEX  OF  THE  CLOSEST  ESTIMAT  (NONPARA  PROCEDURE) 

DELTA  -  NONCENTRALITY  OF  T  DISTRIBUTION 

DF  -  DEGREE  OF  FREEDOM  OF  T  DISTRIBUTION 

EPS  -  WIDTH  OF  THE  CONFIDENCE  BOUND 

EROR  -  EROR  OF  ESTIMATION  (APPROXI  PROCEDURE) 

ERORN  -  EROR  OF  ESTIMATION  (NONPARA  PROCEDURE) 

ERBAR  -  MEAN  OF  EROR  (APPROXI  PROCEDURE) 

ERBARN  -  MEAN  OF  EROR  (NONPARA  PROCEDURE) 

ERSSQ  -  SUM  OF  SQUARE  OF  EROR  (APPROXI  PROCEDURE) 

ERSSQN  -  SUM  OF  SQUARE  OF  EROR  (NONPARA  PROCEDURE) 

ERSUM  -  SUM  OF  EROR  (APPROXI  PROCEDURE) 

ERSUMN  -  SUM  OF  EROR  (NONPARA  PROCEDURE) 

ERSV  -  SAMPLE  VARIANCE  OF  EROR  (APPROXI  PROCEDURE) 

ERSVN  -  SAMPLE  VARIANCE  OF  EROR  (NONPARA  PROCEDURE) 

K  -  STATISTIC  TO  ESTIMATE  DELTA 

LNORM  -  RANDOM  NUMBER  GENERATOR  FOR  NORMAL  VARIATES 

M  -  SAMPLE  SIZE  OF  Y  RANDOM  VARIABLE 

MUX  -  POPULATION  MEAN  OF  X  RANDOM  VARIABLE 

HUY  -  POPULATION  MEAN  OF  Y  RANDOM  VARIABLE 

N  -  SAMPLE  SIZE  OF  X  RANDOM  VARIABLE 

NU  -  THE  SMALLER  OF  SAMPLE  SIZES 

R  -  REAL  RELIABILITY 

RB  -  CONFIDENCE  BOUND  OF  THE  RELIABILITY 

REP  -  REPETITION  OF  SIMULATION 

RLHAT  -  LOWER  CONFIDENCE  LIMIT  OF  RELIABILITY 

RTILD  -  POINT  ESTIMATOR  OF  THE  RELIABILITY 

SIGMAX  -  POPULATION  STANDARD  DEVIATION  OF  X 

SIGMAY  -  POPULATION  STANDARD  DEVIATION  OF  Y 

SP  -  POOLED  SAMPLE  VARIANCE  OF  X  AND  Y 

SUMSQX  -  SUM  OF  SQUARES  OF  X 


40 


* 


SUMSQY 

-  SUM  OF  SQUARES  OF  Y 

* 

SUMX 

-  SUM  OF  X 

* 

SUMY 

-  SUM  OF  Y 

* 

TIN 

-  IMSL  FUNCTION  TO  COMPUTE  PERCENTILE  OF  T  DIST. 

* 

XBAR 

-  AVERAGE  OF  X 

YBAR 

-  AVERAGE  OF  Y 

* 

* 


REQUIRED  EXTERNAL  FUNCTIONS:  LNORM,  TIN,  ANORDF,  ANORIN 


. “'f »'» V%  V* V» »'» V* #V  «*»  **f  i1?  «*»  »W»  V*  •*»•** V*  V»  »'»#'? Vf Vc Vr Vf V» Vr Vf V*  V? 'jV Vr **f *'*  #*»  #'* 1 


*>v 

Vf 


•***;V**i: 


INTEGER  REP,  CASE 
REAL  ALPHA,  TWO 
PARAMETER(ALPHA=0.  05) 

PARAMETER(TW0=2. 0) 

PARAHETER(REP=1000) 

PARAMETER(CASE=18) 

INTEGER  I,  II,  J,  U,  V,  XSEED,  YSEED,  N,  M,  CLOSE, 

+  A,  B,  Al,  Bl,  CLOSEN 

REAL  MUX,  MUY,  XBAR,  YBAR,  SIGMAX,  SIGMAY,  DF,  R,  RLHAT(REP) , 

+  K,  DELTA,  X( 120) ,  Y(120),  Xl(120),  Yl(120),  SUMSQX, 

+  SUMSQY,  SUMX ,  SUMY,  TIN,  ANORDF,  RN,  RM,  DIFF,  SIGMA,  TEMP, 
+  ER,  ERSUM,  ERSSQ,  ERBAR,  ERSV, 

+  BIGU,  RTILD,  NU,  EPS,  ANORIN, ERN,  ERSUMN,  ERSSQN,  ERBARN, 

+  ERSVN,  DFF,  RB(REP) 


CALL  EXCMSC'FILEDEF  12  DISK  SEQ4  DATA  Al’) 

CALL  EXCMSCFILEDEF  18  DISK  AEQ4  DATA  Al') 

DO  2000  1=1,  CASE 

READ  (12,100)  XSEED,  YSEED,  MUX,  MUY,  SIGMAX,  SIGMAY,  N,  M,  R 
100  FORMAT  (15, IX, 15, IX, F5.  1, IX, F7.  3, IX, F4.  1,1X,F4.  1 , IX, 12, IX, 12, IX, 
+F5. 3) 

DF  =  REAL(N+M-2) 

RN  =  REAL(N) 

RM  =  REAL(M) 

ERSUMN  =  0.  0 
ERSSQN  =  0.  0 
ERSUM  =0.0 
ERSSQ  =  0.  0 
DO  1000  J=l,  REP 

CALL  LNORM( XSEED,  X,  N,  2,  0) 

CALL  LNORM( YSEED,  Y,  M,  2,  0) 

SUMSQX  =  0.  0 
SUMSQY  =  0.  0 
SUMX  =  0.  0 
SUMY  =  0.  0 
DO  200  U=  1,  N 

X1(U)  =  X(.U)  *  SIGMAX  +  MUX 
SUMSQX=  SUMSQX  +  X1(U)  *  X1(U) 

SUMX  =  SUMX  +  X1(U) 

200  CONTINUE 

XBAR=  SUMX  /  RN 
DO  300  V=l,  M 

Y1(V)  =  Y(V)  *  SIGMAY  +  MUY 
SUMSQY  =  SUMSQY  +  Y1(V)  *  Y1(V) 
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SUMY  =  SUMY  +  Y1(V) 
300  CONTINUE 


*  PROCEDURE  FOR  PARAMETRIC 
YBAR  =  SUMY  /  RM 

SIGMA  =  SQRT(  (SUMSQX  -  RN*XBAR*XBAR  +  SUMSQY  -  RM*YBAR*YBAR) 
+  /  DF) 

K  =  MAX(  (XBAR  -  YBAR)  /  SIGMA,  0.0  ) 

DELTA  =  K  -  SQRT(  ( RN+RM ) / ( RN*RM )  +  K*K  /  (2. 0*(RN+RM-2. 0))) 
+  *  TIN( 1.0 -ALPHA, DF) 

RLHAT(J)  =  ANORDFC DELTA/ SQRT( TOO)) 

ER  =  R  -  RLHAT(J) 

ERSUM  =  ERSUM  +  ER 
ERSSQ  =  ERSSQ  +  ER  *  ER 


*  PROCEDURE  FOR  NONPARAMETRIC 
BIGU  =  0. 0 
DO  500  A  =  1,  N 
DO  400  B  =  1,  M 

IF  (X1(A)  .GT.  Y1(B))  BIGU  =  BIGU  +1.0 
400  CONTINUE 

500  CONTINUE 

RTILD  =  BIGU  /  (RN  *  RM) 

NU  =  MIN  (RN,  RM) 

EPS  =  1.  0  /  SQRT(4.  0  *  NU)  *  AN0RIN(1.  0  -  ALPHA) 
RB( J)  =  RTILD  -  EPS 
ERN  =  R  -  RB( J) 

ERSUMN  =  ERSUMN  +  ERN 
ERSSQN  =  ERSSQN  +  ERN  *  ERN 
1000  CONTINUE 


ERBARN  =  ERSUMN  /  REAL(REP) 

ERSVN  =  (  ERSSQN  -  REAL(REP)  *  ERBARN  *  ERBARN  )  /  REAL(REP-l) 
DFF  =2.0 

DO  1300  A1  =  1,  REP 

DO  1200  B1  =  A1  +  1,  REP 

IF  (RB(B1)  .LT.  RB(Al))  THEN 
TEMP  =  RB(A1) 

RB(A1)  =  RB(B1) 

RB(B1)  =  TEMP 
ENDIF 

1200  CONTINUE 

IF  (((R-RB(Al))  .GE.  0. IE-6)  .AND. 

+  ((R-RB(Al))  .LE.  DFF))  THEN 

DFF  =  R  -  RB(A1) 

CLOSEN  =  A1 
ENDIF 
1300  CONTINUE 


*  PROCEDURE  OF  PARAMETRIC 

ERBAR  =  ERSUM  /  REAL(REP) 

ERSV  =  (  ERSSQ  -  REAL(REP)  *  ERBAR  *  ERBAR  )  /  REAL(REP-l) 
DIFF  =  2.  0 
DO  1800  11=1,  REP 

DO  1500  J=I1+1,  REP 

IF  (RLHAT(J)  .LT.  RLHAT(Il))  THEN 


TEMP  =  RLHAT(IT) 

RLHAT(Il)  =  RLHAT(J) 

RLKAT(J)  =  TEMP 
ENDIF 

1500  CONTINUE 

IF  (((R-RLHAT(Il))  .  GE.  0. IE-6)  .AND. 

+  ((R-RLHAT(Il))  .LE.  DIFF))  THEN 
DIFF  =  R  -  RLHAT(Il) 

CLOSE  =11 
ENDIF 
1800  CONTINUE 

WRITE  (18,1900)  I,  N,  M,  R,  SIGMAX 

WRITE  (18,1910)  RLHAT(NINT(REAL(REP)*(1- O-ALPHA))),  ERBAR 
WRITE  (18,1920)  RB(NINT(REAL(REP)*(1.  O-ALPHA))),  ERBARN 
1900  FORMAT( 'CASE:  1 ,12,/,'N:  ',  12,  T16,  'M:  ',  12,  /,  'R:  ',  F4.  3, 

+  T16,  'SIGMA  :  '  F4. 1) 

1910  FORMAT('<  PARAMETRIC  >f,  /,  ’RLHAT:  ',  F5.4,  T16,  'ERROR:  ',  F5. 4) 
1920  FORMAT( ' <  NONPARAMETRIC  >',  /,  'RLHAT:  ',  F5.4,  T16,  'ERROR:  ', 

+  F5.4,//) 

2000  CONTINUE 
STOP 
END 
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APPENDIX  D.  FORTRAN  CODE  FOR  INTERVAL  ESTIMATION 

PROCEDURE  -  NORMAL  UNEQUAL  VAR  IANCES  CASE 

PROGRAM  UNEQSM 


'iVVf^VVrxVfiV^iVVciVVrVcicib'cVrj'f^VfiVVf^nV'fyc 


THIS  PROGRAM  IS  TO  VALIDATE  THE  LOWER  CONFIDENCE  LIMIT  APPROXI-  * 
MATION  PROCEDURE  FOR  P(  X  >  Y  ),  WHERE  X,  Y  ARE  INDEPENDENTLY  * 
NORMALLY  DISTRIBUTED  WITH  UNKNOWN  MEANS  AND  UNKNOWN  AND  UNEQUAL  * 
VARIANCES 


ALPHA 

ANORDF 

ASVX 

ASVY 

CASE 

ERROR 

EKRBAR 

ERRSSQ 

ERRSUM 

ERRSV 

CLOSE 

DELTAH 

DF 

LNORM 

M 

MUX 

MUY 

N 

R 

REP 

RLHAT 

SDHAT 

SIGMAX 

SIGMAY 

SUMSQX 

SUHSQY 

SUMX 

SUMY 

SVX 

SVY 

T 

TIN 

VARHAT 

XBAR 

YBAR 


NOMINAL  CONFIDENCE  LEVEL 

IMSL  FUNCTION  FOR  NORMAL  PROBABILITY  * 

SAMPLE  VARIANCE  OF  X  DEVIDED  BY  SAMPLE  SIZE 

SAMPLE  VARIANCE  OF  Y  DEVIDED  BY  SAMPLE  SIZE  * 

NUMBER  OF  TESTING 

ERROR  OF  ESTIMATION 

MEAN  OF  ERROR 

SUM  OF  SQUARE  OF  ERROR 

SUM  OF  ERROR 

SAMPLE  VARIANCE  OF  ERROR  * 

INDEX  OF  THE  CLOSEST  ESTIMAT  * 

DELTAHAT;  A  ESTIMATOR  OF  DELTA 

DEGREES  OF  FREEDOM  OF  DELTAHAT 

RANDOM  NUMBER  GENERATOR  FOR  NORMAL  VARIATES 

SAMPLE  SIZE  OF  RANDOM  VARIABLE  X  * 

POPULATION  MEAN  OF  X 

POPULATION  MEAN  OF  Y 

SAMPLE  SIZE  OF  RANDOM  VARIABLE  Y  * 

REAL  RELIABILITY 

REPETITION  OF  SIMULATIONS 

LOWER  CONFIDENCE  LIMIT  OF  RELIABILITY 

SAMPLE  STANDARD  DEVIATION  OF  THE  DELTAHAT 

STANDARD  DEVIATION  OF  X  * 

STANDARD  DEVIATION  OF  Y  * 

SUM  OF  SQUARE  OF  RANDOM  SAMPLE  OF  X 

SUM  OF  SQUARE  OF  RANDOM  SAMPLE  OF  Y 

SUM  OF  RANDOM  SAMPLE  OF  X  * 

SUM  OF  RANDOM  SAMPLE  OF  Y 

SAMPLE  VARIANCE  OF  X  * 

SAMPLE  VARIANCE  OF  Y 

PERCENTILE  OF  THE  T  DISTRIBUTION  * 

IMSL  FUNCTION  TO  COMPUTE  PERCENTILE  OF  T  DISTRIBUTION  * 
SAMPLE  VARIANCE  OF  DELTAHAT 
SAMPLE  MEAN  OF  X 

SAMPLE  MEAN  OF  Y  * 


REQUIRED  EXTERNAL  FUNCTIONS:  LNORM,  TIN,  ANORDF 


INTEGER  REP,  CASE 
REAL  ALPHA 
PARAMETER(REP=1000) 

PAR AMETER( CASE=1 8 ) 

PARAMETER(ALPHA=0. 05) 

INTEGER  I,  11,  J,  U,  V,  XSEED.,  YSEED,  N,  M,  CLOSE 
REAL  RM,  RN,  MUX,  MUY,  XBAR,  YBAR,  SIGMAX,  SIGMAY,  R, 

+  X(100),  Y(100) ,  XI,  Yl,  SUMSQX,  SUMSQY,  SUMX,  SUMY, 

+  SVX,  SVY,  ASVX,  ASVY,  DF1,  DF,  T,  TIN,  DELTAH,  VARHAT, 

4-  SDHAT,  RLHAT(REP) ,  ANORDF,  TEMP,  DIFF,  ERROR,  ERRSUM,  ERRSSQ, 

4-  ERRBAR,  ERRSV 


CALL  EXCMSC 'FILEDEF  12  DISK  SETUNQ  DATA  Al') 

CALL  EXCMS( 'FILEDEF  18  DISK  0PT2  DATA  Al') 

WRITE  (18,*)  'UNEQUAL  VARIANCES' 

DO  2000  1=1,  CASE 

READ  (12,100)  XSEED,  YSEED,  MUX,  MUY,  SIGMAX,  SIGMAY,  N,  M,  R 
100  FORMAT  (15, IX, 15, IX, F5.  1,1X,F7. 3,1X,F4.  1,1X,F4.  1, IX, 12, IX, 12, IX, 
+F5. 3) 

RN  =  REAL(N) 

RM  =  REAL(M) 

ERRSUM  =  0.  0 
ERRSSQ  =  0.  0 

*  <  TRANSFORMATION  OF  X,  Y  TO  DESIRED  PROPERTIES  > 

DO  1000  J=l,  REP 

CALL  LNORM( XSEED,  X,  N,  2,  0) 

CALL  LNORMC YSEED,  Y,  M,  2,  0) 

SUMSQX  =  0.  0 
SUMSQY  =  0.  0 
SUMX  =  0.0 
SUMY  =  0.0 
DO  200  U=  1,  N 

XI  =  X(U)  *  SIGMAX  4-  MUX 
SUMSQX=  SUMSQX  4-  XI  *  XI 
SUMX  =  SUMX  4-  XI 
200  CONTINUE 

XBAR=  SUMX  /  RN 
DO  300  V=l,  M 

Yl  =  Y(V)  *  SIGMAY  4-  MUY 
SUMSQY  =  SUMSQY  4-  Yl  *  Yl 
SUMY  =  SUMY  4-  Yl 
300  CONTINUE 

*  <  COMPUTE  CONFIDENCE  LIMIT  OF  RELIABILITY  > 

YBAR  =  SUMY  /  RM 

SVX  =  (SUMSQX  -  RN  *  XBAR  *  XBAR)  /  (RN  -  1.0) 

SVY  =  (SUMSQY  -  RM  *  YBAR  *  YBAR)  /  (RM  -  1.0) 

ASVX  =  SVX  /  RN 
ASVY  =  SVY  /  RM 

DF1  =  (SVX  4-  SVY)  *  (SVX  4-  SVY)  /  (  SVX*SVX  /  (RN-1.0)  4- 
4-  SVY*SVY  /  (RM-1.0)  ) 

DF  =  ANINT(DFl) 

T  =  TIN(  1. O-ALPHA,  DF) 

DELTAH  =  (XBAR  -  YBAR)  /  SQRT(SVX  +  SVY) 
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+ 

+ 

+ 


VARHAT  =  (ASVX  +  ASVY)  /  (SVX  +  SVY)  + 

(XBAR  -  YBAR)  *  (XBAR  -  YBAR)  / 

(  2.0* (SVX  +  SVY)**3  )  * 

(  ( SVX**2 )  /  (RN-1. 0)  +  (SVY**2)  /  (RM-1.0)  ) 

SDHAT  =  SQRT( VARHAT) 

RLHAT(J)  =  AN0RDF(DELTAH  -  T  *  SDHAT) 

<  COMPUTE  THE  MEAN  AND  VARIANCE  OF  ESTIMA1.  ‘N  ERROR  > 

ERROR=  R  -  RLHAT(J) 

ERRSUM  =  ERRSUM  +  ERROR 
ERRSSQ  =  ERRSSQ  +  ERROR  *  ERROR 
ERRBAR  =  ERRSUM  /  REAL(REP) 

ERRSV  =  (  ERRSSQ  -  REAL(REP)  *  ERRBAR  *  ERRBAR  )  /  REAL(REP-l) 
1000  CONTINUE 

*  <  SORT  CONFIDENCE  LIMITS  IN  ASCENDING  ORDER  > 


DIFF  =  2.  0 
DO  1800  11=1,  REP 

DO  1500  J=I1+1,  REP 

IF  (RLHAT(J)  .LT.  RLHAT(Il))  THEN 
TEMP  =  RLHAT(Il) 

RLHAT(Il)  =  RLHAT(J) 

RLHAT(J)  =  TEMP 
END  IF 

1500  CONTINUE 

Vc 

*  <  FIND  THE  CLOSEST  CONFIDENCE  LIMIT  ESTIMATE  > 


1800 


1900 


2000 


IF  (((R-RLHAT(Il))  . GE.  0. IE-6)  .AND. 

+  ((R-RLHAT(Il))  .LE.  DIFF))  THEN 
DIFF  =  R  -  RLHAT(Il) 

CLOSE  =11 
ENDIF 
CONTINUE 

WRITE  (18,1900)  I,  MUX,  N,  SIGMAX,  MUY,  M,  SIGMAY,  R, 

+  RLHAT(NINT(REAL(REP)*( 1. 0 "ALPHA) ) ) ,  RLHAT( CLOSE) , 

+  REAL( CLOSE)/ 1000.  0,  ERRBAR,  ERRSV 

FORMATC SIMULATION:  ’,12,/,' MUX:  ' ,F5. 1 ,T16 , ' N:  ',12, 

+  T35,' SIGMAX:  ' ,F4. 1,/, 'MUY:  1 ,F7.  3,T16, 'M:  ',12, 

+  T35, 'SIGMAY:  ' ,F4. 1,/  'TRUE  R  :  ',F7.5, 

+  T35, 'RLHAT:  ' ,F7. 5 ,/, 'CLOSEST  RLHAT:  ',F7.5, 

+  T35 , ' TRUE  CONFIDENCE  LEVEL:  ' ,F5 . 3 , / , 

+  'MEAN  ERROR:  ',F7.5,/, 

+  'VARIANCE  OF  ERROR:  ’,F7.5,///) 

CONTINUE 

STOP 

END 


46 


APPENDIX  E.  SIMULATION  PARAMETER  SETS  FOR  UNEQUAL 


VARIANCES  CASE 


FILE: 

SETUNQ  DATA 

(FOR  APPROXIMATE  PROCEDURE) 

16807 

93943 

10.0 

7.  133 

1.0 

2.  0 

10 

20 

.900 

1 

16807 

93943 

10.0 

7.  133 

1.0 

2.0 

25 

35 

.  900 

2 

16807 

93943 

10.  0 

7.  133 

1.  0 

2.0 

75 

50 

.  900 

3 

16807 

93943 

300.  0 

247. 142 

10.  0 

40.  0 

10 

20 

.  900 

4 

16807 

93943 

300.  0 

247. 142 

10.0 

40.0 

25 

35 

.  900 

5 

16807 

93943 

300.  0 

247. 142 

10.  0 

40.  0 

75 

50 

.900 

6 

16807 

93943 

10.  0 

6.-322 

1.  0 

2.  0 

10 

20 

.  950 

7 

16807 

93943 

10.  0 

6.  322 

1.  0 

2.0 

25 

35 

.  950 

8 

16807 

93943 

10.-0 

6.  322 

1.0 

2.0 

75 

50 

.950 

9 

16807 

93943 

300.  0 

232. 175 

10.  0 

40.  0 

10 

20 

.  950 

10 

16807 

93943 

300.  0 

232. 175 

10.  0 

40.  0 

25 

35 

.  950 

11 

16807 

93943 

300.  0 

232. 175 

10.  0 

40.  0 

75 

50 

.  950 

12 

16807 

93943 

10.  0 

4.  799 

1.  0 

2.0 

10 

20 

.  990 

13 

16807 

93943 

10.  0 

4.  799 

1.  0 

2.  0 

25 

35 

.  990 

14 

16807 

93943 

10.  0 

4.  799 

1.  0 

2.  0 

75 

50 

.  990 

15 

16807 

93943 

300.  0 

204. 097 

10.0 

40.  0 

10 

20 

.990 

16 

16807 

93943 

300.0 

204. 097 

10.0 

40.  0 

25 

35 

.  990 

17 

16807 

93943 

300.  0 

204. 097 

10.0 

40.0 

75 

50 

.990 

18 

X 

Y 

X 

Y 

X 

Y 

SEED 

SEED 

MEAN1 

MEAN 

STDV 

STDV 

N 

M 

R 

* 
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FILE: 

SUQ7  DATA  (FOR  COMPARISON  SIMULATIONS). 

16807 

93943 

10.0 

7.  133 

1.  0 

2.  0 

10 

15 

.  900 

1 

16807 

93943 

10.0 

7.  133 

1.  0 

2.-0 

70 

35 

.900 

2 

16807 

93943 

10.0 

7.  133 

1.  0 

2.  0 

90 

90 

.  900 

3 

16807 

93943 

300.0 

247.  142 

10.  0 

40.0 

10 

15 

.  900 

4 

16807 

93943 

300.  0 

247.  142 

10.  0 

40.0 

70 

35 

.  900 

5 

16807 

93943 

300.  0 

247. 142 

10.  0 

40.0 

90 

90 

.  900 

6 

16807 

93943 

10.0 

6.322 

1.  0 

2.  0 

10 

15 

.950 

7 

16807 

93943 

10.  0 

6.322 

1.  0 

2.  0 

70 

35 

.  950 

8 

16807 

93943 

10.0 

6.322 

1.  0 

2.  0 

90 

90 

.  950 

9 

16807 

93943 

300.  0 

232.  175 

10.0 

40.  0 

10 

15 

.  950 

10 

16807 

93943 

300.0 

232.  175 

10.  0 

40.  0 

70 

35 

.950 

11 

16807 

93943 

300.0 

232.  175 

10.0 

40.0 

90 

90 

.950 

12 

16807 

93943 

10.0 

4.  799 

1.  0 

2.0 

10 

15 

.990 

13 

16807 

93943 

10.0 

4.  799 

1.  0 

2.-0 

70 

35 

.990 

14 

16807 

93943 

10.  0 

4.  799 

1.  0 

2.  0 

90 

90 

.  990 

15 

16807 

93943 

300.  0 

204. 097 

10.  0 

40.  0 

10 

15 

.  990 

16 

16807 

93943 

300.  0 

204. 097 

10.  0 

40.  0 

70 

35 

.  990 

17 

16807 

93943 

300.0 

204. 097 

10.  0 

40.0 

90 

90 

.  990 

18 

IS 

15 

F5.  1 

F7.  3 

F4.  1 

F4.  1 

12 

12 

F5.  3 

X 

Y 

X 

Y 

X 

Y 

N 

M 

R 

SEED 

SEED 

MEAN 

MEAN 

STDV 

STDV 

4S 


APPENDIX  F.  FORTRAN  CODE  FOR  SIMULTANEOUS  COMPARISON 
SIMULATION  OF  APPROXIMATE  PROCEDURE  VS. 


NONPARAMETRIC  PROCEDURE  -  UNEQUAL  VARIANCES 


PROGRAM  C0MUQ7 

* 


***  *  **  ****  *  **  **  *****  ***  *  *  ******  *  **  *  *  **  *  **  *  **  **  **  ****  *****  **  ***  **  ******* 


THIS  PROGRAM  IS  TO  COMPUTE  THE  LOWER  CONFIDENCE  LIMITS  OF 
R  =  P  (X  >  Y)  SIMULTANEOUSLY  FOR  BOTH  THE  APPROXIMATE  PRO-  * 

AND  THE  NONPARAMETRIC  PROCEDURE,  WHERE  X,  Y  ARE  INDEPENDENTLY  * 

NORMALLY  DISTRIBUTED  WITH  UNKNOWN  MEANS  AND  UNKNOWN  AND  UNEQUAL  * 
VARIANCES  * 


* 

ALPHA 

- 

NOMINAL  CONFIDENCE  LEVEL 

* 

ANORDF 

- 

IMSL  FUNCTION  FOR  NORMAL  PROBABILITY 

* 

ANORIN 

- 

IMSL  FUNCTION  FOR  INVERSE  NORMAL  CDF 

* 

ASVX 

- 

SAMPLE  VARIANCE  OF  X  DEVIDED  BY  SAMPLE  SIZE 

* 

ASVY 

- 

SAMPLE  VARIANCE  OF  Y  DEVIDED  BY  SAMPLE  SIZE 

* 

BIGU 

- 

MANN-WHITNEY  STATISTIC 

* 

CASE 

- 

NUMBER  OF  TESTING 

* 

EPS 

- 

WIDTH  OF  THE  CONFIDENCE  BOUND 

* 

EROR 

- 

EROR  OF  ESTIMATION  (APPROXI  PROCEDURE) 

* 

ERORN 

- 

EROR  OF  ESTIMATION  (NONPARA  PROCEDURE) 

* 

ERBAR 

- 

MEAN  OF  EROR  (APPROXI  PROCEDURE) 

* 

ERBARN 

- 

MEAN  OF  EROR  (NONPARA  PROCEDURE) 

* 

ERSSQ 

- 

SUM  OF  SQUARE  OF  EROR  (APPROXI  PROCEDURE) 

* 

ERSSQN 

- 

SUM  OF  SQUARE  OF  EROR  (NONPARA  PROCEDURE) 

* 

ERSUM 

- 

SUM  OF  EROR  (APPROXI  PROCEDURE) 

* 

ERSUMN 

- 

SUM  OF  EROR  (NONPARA  PROCEDURE) 

* 

ERSV 

- 

SAMPLE  VARIANCE  OF  EROR  (APPROXI  PROCEDURE) 

* 

ERSVN 

- 

SAMPLE  VARIANCE  OF  EROR  (NONPARA  PROCEDURE) 

* 

CLOSE 

- 

INDEX  OF  THE  CLOSEST  ESTIMAT  (APPROXI  PROCEDURE) 

* 

CLOSEN 

- 

INDEX  OF  THE  CLOSEST  ESTIMAT  (NONPARA  PROCEDURE) 

* 

DELTAH 

- 

DELTAHAT;  A  ESTIMATOR  OF  DELTA 

* 

DF 

- 

DEGREES  OF  FREEDOM  OF  DELTAHAT 

* 

LNORM 

- 

RANDOM  NUMBER  GENERATOR  FOR  NORMAL  VARIATES 

* 

M 

- 

SAMPLE  SIZE  OF  RANDOM  VARIABLE  X 

* 

MUX 

- 

POPULATION  MEAN  OF  X 

* 

MUY 

- 

POPULATION  MEAN  OF  Y 

* 

N 

- 

SAMPLE  SIZE  OF  RANDOM  VARIABLE  Y 

* 

NU 

- 

THE  SMALLER  OF  SAMPLE  SIZES 

* 

R 

- 

REAL  RELIABILITY 

* 

RB 

- 

CONFIDENCE  BOUND  OF  THE  RELIABILITY 

* 

REP 

- 

REPETITION  OF  SIMULATIONS 

* 

RLHAT 

- 

LOWER  CINFIDENCE  LIMIT  OF  RELIABILITY 

* 

RTILD 

- 

POINT  ESTIMATOR  OF  THE  RELIABILITY 

* 

SDHAT 

- 

SAMPLE  STANDARD  DEVIATION  OF  THE  DELTAHAT 

* 

SIGMAX 

- 

STANDARD  DEVIATION  OF  X 

* 

SIC-MAY 

- 

STANDARD  DEVIATION  OF  Y 

* 

SUMSQX 

- 

SUM  OF  SQUARE  OF  RANDOM  SAMPLE  OF  X 

* 

* 

* 

* 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


* 


* 

* 

* 

* 

* 

* 


* 

* 


* 

* 

* 

* 


* 

* 

* 

* 


* 

* 
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* 

SUMSQY 

- 

SUM  OF  SQUARE  OF  RANDOM  SAMPLE  OF  Y 

4c 

SUMX 

- 

SUM  OF  RANDOM  SAMPLE  OF  X 

4c 

SUMY 

- 

SUM  OF  RANDOM  SAMPLE  OF  Y 

4c 

SVX 

- 

SAMPLE  VARIANCE  OF  X 

* 

SVY 

- 

SAMPLE  VARIANCE  OF  Y 

4c 

T 

- 

PERCENTILE  OF  THE  T  DISTRIBUTION 

4c 

TIN 

- 

IMSL  FUNCTION  TO  COMPUTE  PERCENTILE  OF  T  DISTRIBUTION 

4c 

VARHAT 

- 

SAMPLE  VARIANCE  OF  DELTAHAT 

4c 

XBAR 

- 

SAMPLE  MEAN  OF  X 

4c 

* 

YBAR 

- 

SAMPLE  MEAN  OF  Y 

* 

* 

REQUIRED  EXTERNAL  FUNCTIONS:  LNORM,  TIN,  ANORDF,  ANORIN 

**iV5Wc*****5V*5ViViV*5V*3V5V***5Wf5V3Wr*5V*5Wf!V*>WrjWcs'c5WciV-'r************5'r5WfVr5WrVr*******7W: 


INTEGER  REP,  CASE 

REAL  ALPHA 

PAR AMETER ( RE  P= 1000) 

PARAMETER(CASE=18) 

PARAMETER( ALPHA=0.  05) 

INTEGER  I,  II,  J,  U,  V,  XSEED,  YSEED,  N,  M,  CLOSE, 

+  A,  B,  Al,  Bl,  CLOSEN 

REAL  RM,  RN,  MUX,  MUY,  XBAR,  YBAR,  SIGMAX,  SIGMAY,  R, 

+  X(100),  Y(100),  Xl(lOO),  Yl(lOO),  SUMSQX,  SUMSQY,  SUMX,  SUMY, 

+  SVX,  SVY,  ASVX,  ASVY,  DF1,  DF,  T,  TIN,  DELTAH,  VARHAT, 

+  SDHAT,  RLHAT(REP) ,  ANORDF,  TEMP,  DIFF,  EROR,  ERSUM,  ERSSQ, 

+  ERBAR,  ERSV, 

+  BIGU,  RTILD,  NU,  EPS,  ANORIN,  ERN,  ERSUMN,  ERSSQN,  ERBARN, 

+  ERSVN ,  OFF,  RB(REP) 


CALL  EXCMS( 'FILEDEF  12  DISK  SUQ7  DATA  Al') 

CALL  EXCMS( 'FILEDEF  18  DISK  AUQ7  DATA  Al') 

WRITE  (IS,*)  'UNEQUAL  VARIANCES' 

DO  2000  1=1,  CASE 

READ  (12,100)  XSEED,  YSEED,  MUX,  MUY,  SIGMAX,  SIGMAY,  N,  M,  R 
100  FORMAT  (15, IX, 15, IX, F5. 1,1X,F7.  3,1X,F4. 1,1X,F4. 1, IX, 12, IX, 12, IX, 
+F5. 3) 

RN  =  REAL(N) 

RM  =  REAL(M) 

ERSUM  =  0.  0 
ERSSQ  =  0.  0 
ERSUMN  =0.0 
ERSSQN  =  0.  0 
DO  1000  J=l,  REP 

CALL  LNORM( XSEED,  X,  N,  2,  0) 

CALL  LNORM( YSEED,  Y,  M,  2,  0) 

SUMSQX  =  0.0 
SUMSQY  =  0.  0 
SUMX  =  0. 0 
SUMY  =  0. 0 
DO  200  U=  1,  N 

X1(U)  =  X(U)  *  SIGMAX  +  MUX 
SUMSQX=  SUMSQX  +  X1(U)  *  X1(U) 

SUMX  =  SUMX  +  X1(U) 

200  CONTINUE 

XBAR=  SUMX  /  RN 


DO  300  V=1 ,  M 

Y1(V)  =  Y(V)  *  SIGMAY  +  MUY 
SUMSQY  =  SUMSQY  +  Y1(V)  *  Y1(V) 

SUMY  =  SUMY  +  Yl( V) 

300  CONTINUE 

YBAR  =  SUMY  /  RM 

SVX  =  (SUMSQX  -  RN  *  XBAR  *  XBAR)  /  (RN  -  1.0) 

SVY  =  (SUMSQY  -  RM  *  YBAR  *  YBAR)  /  (RM  -  1.0) 

ASVX  =  SVX  /  RN 
ASVY  =  SVY  /  RM 

DF1  =  (SVX  +  SVY)  *  (SVX  +  SVY)  /  (  SVX*SVX  /  (RN-1.0)  + 
+  SVY*SVY  /  (RM-1.  0)  ) 

DF  =  ANINT(DFl) 

T  =  TIN(  1.  O-ALPHA,  DF) 

DELTAH  =  (XBAR  -  YBAR)  /  SQRT(SVX  +  SVY) 

VARHAT  =  (ASVX  +  ASVY)  /  (SVX  4-  SVY)  + 

+  (XBAR  -  YBAR)  *  (XBAR  -  YBAR)  / 

+  (  2.  0*(SVX  +  SVY) **3  )  * 

+  (  (SVX**2)  /  (RN-1.0)  +  (SVY**2)  /  (RM-1.  0)  ) 

SDHAT  =  SQRT( VARHAT) 

RLHAT(J)  =  ANORDF( DELTAH  -  T  *  SDHAT) 

EROR  =  R  -  RLHAT(J) 

ERSUM  =  ERSUM  +  EROR 
ERSSQ  =  ERSSQ  +  EROR  *  EROR 


*  PROCEDURE  FOR  NONPARAMETRIC 
BIGU  =  0. 0 
DO  500  A  =  1,  N 
DO  400  B  =  1,  M 

IF  (X1(A)  .GT.  Y1(B))  BIGU  =  BIGU  +1.0 
400  CONTINUE 

500  CONTINUE 

RTILD  =  BIGU  /  (RN  *  RM) 

NU  =  MIN  (RN,  RM) 

EPS  =  1.  0  /  SQRT(4. 0  *  NU)  *  AN0RIN(1.  0  -  ALPHA) 
RB( J)  =  RTILD  -  EPS 
ERN  =  R  -  RB( J) 

ERSUMN  =  ERSUMN  +  ERN 
ERSSQN  =  ERSSQN  +  ERN  *  ERN 
1000  CONTINUE 


ERBARN  =  ERSUMN  /  REAL(REP) 

ERSVN  =  (  ERSSQN  -  REAL(REP)  *  ERBARN  *  ERBARN  )  /  REAL(REP-l) 
DFF  =  2.  0 

DO  1300  Al  =  1,  REP 

DO  1200  B1  =  Al  +  1,  REP 

IF  (RB(B1)  .LT.  RB(A1))  THEN 
TEMP  =  RB(A1) 

RB(A1)  =  RB(B1) 

RB(31)  =  TEMP 
ENDIF 

1200  CONTINUE 

IF  (((R-RB(Al))  .GE.  0. IE-6)  .AND. 

+  ((R-RB(Al))  .LE.  DFF))  THEN 

DFF  =  R  -  RB(Al) 

CLOSEN  =  Al 
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ENDIF 
1300  CONTINUE 

Vr -  --  --  -- -  --  --  --  -- - - - 

*  PROCEDURE  OF  PARAMETRIC 

ERBAR  =  ERSUM  /  REAL(REP) 

ERSV  =  (  ERSSQ  -  REAL(REP)  *  ERBAR  *  ERBAR  )  /  REAL(REP-l) 

DIFF  =  2.  0 
DO  1800  11=1,  REP 

DO  1500  J=I1+1,  REP 

IF  (RLHAT(J)  .LT.  RLHAT(Il))  THEN 
TEMP  =  RLHAT(Il) 

RLHAT(Il)  =  RLKAT(J) 

RLHAT(J)  =  TEMP 
ENDIF 

1500  CONTINUE 

IF  (((R-RLHAT(Il))  .  GE.  0. IE-6)  .AND. 

+  ((R-RLHAT(ri))  .LE.  DIFF))  THEN 
DIFF  =  R  -  RLHAT(Il) 

CLOSE  =  II 
ENDIF 
1800  CONTINUE 

WRITE  (18,1500)  I,  N,  M,  R,  SIGMAX,  SIGMAY 
WRITE  (18,1910)  RLHAT(NINT(REAL(REP)*( 1.0-ALPHA))),  ERBAR 
WRITE  (18,1920)  RB(NINT(REAL(REP)*( 1.  0-ALPHA) )) ,  ERBARN 
1900  FORMAT( ' CASE:  ',12,/,'N:  12,  T16,  'M:  ’,  12.  /,  'R:  ',  F4.3, 

-f  T16,  'SIGMAX:  1  ,  F4.  T,  T35,  'SIGMAY:  f,F4.  1) 

1910  FORMAT( ' <  PARAMETRIC  >' ,  /,  'RLHAT:  ',  F5.4,  T16,  'ERROR:  ',  F5.4) 
•1920  FORMAT( ' <  NONPARAMETRIC  >' ,  /,  'RLHAT:  ',  F5.4,  T16,  'ERROR:  ', 

+  F5.4,//) 

2000  CONTINUE 
STOP 
END 
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